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I. INTRODUCTION 



There are a number of methods in use today that numeri- 
cally solve for the currents on simple wire antenna struc- 
tures. The majority of these methods utilize an integral 
equation solution to Maxwell's Equations which also match the 
boundary conditions on the structure and satisfy the radiation 
conditions for the radiated fields. One of the most widely 
used implementations of the integral equation method is the 
Numerical Electromagnetics Code (NEC) [Ref 6] . This program 
has been proven to be highly accurate for antenna structures 
that can be modeled using individual wire segments or surface 
patch elements. 

A new method is presented here which does not rely upon 
an integral equation formulation. This technique encloses the 
antenna structure between two parallel ground planes, as an 
approximation to the free space condition, to allow a partial 
field representation using a periodic Fourier series. In the 
initial investigation being reported here, the simple monopole 
radiator is being considered. In such a case, cylindrical 
regions are established around the structure and the harmonic 
expansions of the magnetic and electric fields are matched 
along the regional boundaries to solve for the unknown 
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expansion coefficients. These coefficients are then used to 
compute the current imposed on the surface of the antenna. 

This project began as a search for the viability of 

computing the currents on a top loaded monopole between two 
ground planes using cylindrical harmonic expansions in three 
regions. After a number of weeks of inconclusive results, it 
was decided to look at the simpler case of the monopole 

antenna between two ground planes as presented in Chapter 2 . 

During the validation of the code, as presented in 

Chapter 3 , it was found that the number of terms used in each 
region can be reduced by as much as 90% and still yield 
consistent results for the surface current and input im- 
pedance. However, the representation of the electric field 
along the surface of the antenna is highly inaccurate when the 
number of terms are insufficient to adequately represent the 
field. This phenomenon will be discussed in Chapter 4. 
Additional research revealed that a simpler method could be 
used that would remove the requirement for one complete 

expansion. This is presented in Chapter 5. Finally, con- 
clusions and recommendations are presented in the last 
chapter. In addition, Appendices are used to contain much of 
the detailed descriptions of the computer algorithms and 
resultant programming for this effort. 
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A. METHOD USED 



The general approach is to place the antenna structure 
being analyzed between two perfectly conducting ground planes 
and divide the surrounding area into a number of cylindrical 
regions, using a cylindrical coordinate system, that coincide 
with the natural boundary interfaces of the antenna as shown 
in Figure 1. The upper ground plane is introduced to com- 
pletely close the regions and produces multiple images, each 
identical to the previous image but with opposite polarity. 

Once the regions are defined, we start by obtaining the 
solution to the complex scalar Helmholtz equation. Then, by 





Figure 1 The Unloaded Monopole Antenna Structure 
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using an harmonic expansion of the resulting scalar wave 
function, we write an expression for the vector electric 
and magnetic fields in each region. The tangential component 
of the electric field, as approximated by the truncated 
expansion in Region I, is set equal to the "known” field in 
the gap and set equal to zero along the cylindrical surface 
of the perfectly conducting monopole. The tangential electric 
and magnetic field expansions are then matched across the 
boundary interface between Regions I and II. The resulting 
set of equations is reduced to a system of one equation with 
one unknown by sifting out some of the coefficients using the 
orthogonality principle of Fourier moment integrations. The 
resulting system is then solved for the remaining set of 
unknown coefficients. These coefficients are then used to 
compute the current distribution along the surface of the 
antenna. To validate the source code, the results are then 
compared to those obtained from the Numerical Electromagnetics 
Code (NEC) to ascertain computational accuracy. 

B. ASSUMPTIONS MADE 

In order to simplify the derivations, assumptions were 
made regarding the electric characteristics along the ground 
planes and the surface of the antenna structure, and the 
electric field produced by a driving source of constant 
frequency. First, the ground planes are modeled as perfectly 
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conducting ground planes of infinite dimension. In reality, 
a ground plane can be considered to have infinite dimension, 
for purposes of antenna current calculations, if it is larger 
in dimension than about 10 times the size of the antenna [Ref. 
1]. Additionally, the electric field produced by the driving 
source is considered to be of constant magnitude within the 
gap region. This is a standard source model employed in 
integral equation modeling of antenna structures, [Ref. 1]. 
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II. COMPUTATIONAL PROCEDURE 



This chapter will present the mathematical background 
required to solve for the currents and the input impedance on 
a simple monopole antenna structure using cylindrical harmonic 
expansions. The final solution will be reduced to a system 
of one set of unknowns which will be solved using a digital 
computer system. An overview of the software code used, 
including error detection, hardware requirements and a brief 
discussion of the input and output file structures will also 
be presented. Details of the algorithm implementation are 
relegated to Appendices. 

A. THEORETICAL TECHNIQUE 

The antenna structure shown in Figure 1 (page 3) is 
orientated within the cylindrical coordinate system as shown 
in Figure 2. We begin this development with the complex 
scalar wave equation, known as the Helmholtz equation, given 
in cylindrical coordinates as [Refs. 2, 3] 

2 2 
V y ( p,z ) + k y(p,z ) = 0 

{Eq. 1-1} 

where is used to denote the electric type Hertz potential 
and not the magnetic flux [Ref. 3]. Expanding this equation, 
and noting that the structure and the fields are assumed to 
be axisymmetric, one obtains the following: 
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{Eq. 1-2} 

Utilizing axisymmetry, the transverse magnetic (TM) 
solutions are given in cylindrical coordinates as [Ref. 3] 
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{Eq. 1-3} 



{Eq. 1-4} 
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The following boundary conditionons must be satisfied: 

• The electric field tangential to the ground planes must 
be zero (E„ = 0) . 

• The electric field tangential to the top of the antenna 
structure must be zero. 

• The tangential electric field along the antenna structure 
must be zero (E z = 0) . 

• The electric field along the gap distance d must be a 
constant value (E z = -V Q /d) . 

• The electric and magnetic fields must be continuous 
across the boundary between the two regions above the 
antenna structures (E z , H^) . 

• The radiated fields in Region I must be purely outbound 
at large distances from the antenna and must satisfy the 
radiation condition that E z = rj 0 H^, where rj 0 = 377n. 

Using separation of variables, where the product solution 
has the form \j> = R(p)Z(z), one obtains the following equa- 
tions: [Ref. 3] 



The first equation is recognized as Bessel's equation of order 
n. The solutions to Bessel's equation are represented by any 
of several special functions, including [Ref. 3] 




{Eq. 1-5} 



J„(V)» Y n(V)» H n <1, (M)» H n <2> (V) < E( 5- 1_ 6 > 
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where J n (k p p) and Y n (k p p) represent Bessel functions of the 

first and second kind of order n, and H n a> (/c p p) and H„ (2) (k p p) 

are Hankel functions of the first and second kind of order n, 
(2) 

where H,, = J n ± j Y n . As any two of these functions are 

linearly independent of each other, the solution to Bessel's 
equation can be represented as a linear combination of any two 
of these functions. The solution to the Z equation is a 
linear combinations of harmonic functions of the form e jkz and 
e‘ Jkz . 

Using these results, the desired solutions to the Helmho- 
ltz equation becomes 



where 



Y( A?) = 



J o(k p p) 
H 0 °\k p p) 



j sin (k z z) ) 
\cos (k z z) J 



{Eq. 1-7} 



k 



2 

p 




= k 



2 

o 



2 

= co pe 



(Eq. 1-8} 

and where k 0 is defined as the wave number, co the angular 
frequency, p the permeability of a vacuum ( 4 7r x 1(T 7 H/m) and 
c the permittivity of a vacuum (8.854 x 10’ 12 F/m) . Acceptable 
quantized values of k z are found by substituting the product 
solution into Equation 1-3 and setting E p = 0 on the upper 
ground planes in Region I while enforcing E p = 0 on the upper 
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ground plane and the top of the antenna in Region II. This 
gives, respectively, 



k, = ™ 


in Region I 


{Eq. 1-9. a} 


n l 






K = n7r 


in Region II 


{Eq. 1-9. b} 


n l-h 







It should be noted here that the trigonometric terms used in 
Equation 1-7 will uniquely satisfy the first three boundary 
conditions listed previously for each term in the series. 

For large arguments, the Bessel and Hankel functions 
listed above closely resemble harmonic functions. As seen in 
Table I [Ref* 3], J n (k p p) and Y n (k p p) are analogous to cosine 
and sine functions respectively, while H n tl> ( A: p p ) and H n (2) (k p p) 
are related to complex exponential functions. Therefore, 
Bessel functions can be used to represent standing waves, 
while the Hankel functions describe travelling waves; specifi- 
cally Hn' 15 represents inbound waves while H n (2) represents 
outbound waves. Using these properties, the potential field 
in Region I can be written as the product of Hankel and cosine 
functions where the cosine term implicitly satisfies the 
tangential electric and magnetic field boundary conditions. 

Because the higher spectral bandwidth of the potential 
field generated by the gap voltage was expected to require a 
larger number of terms in the Fourier series expansion than 
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Table I PROPERTIES OF SOLUTIONS TO BESSEL'S EQUATION 
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that of the quasi-singular field at the end of the antenna, 
the total field in Region I was separated into two parts as 
shown below 



¥ (P.Z ) = X ° n H 0 <2) (V n P ) COS (^J + ^ b m H o\ P ) COS ( : 



inn z 



l 



n = 0 



m = 0 



= ¥ a + n 



(Eq. 1-10 } 



where 






{Eq. 1-11} 

In Region II, the potential field is constrained as a 
standing wave and can be expressed as 



V / (P» Z )= X c n J o( ll nP) c °s 



n = 0 



nn 



M) 



{Eq. 1-12} 



where 



u 




{Eq. 1-13} 



and q = l - h. 



Using Equation 1-3, the electric field component at the end 
of the antenna is approximated by the truncated expansion 
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A! = 0 



! 




0 , 



for 0 <z <h 



{Eq. 1-14} 

and the electric field due to the gap voltage is given by 



(Eq. 1-15} 

where N x represents the number of terms to be used to in 
Region I, N 2 for Region II, and for the electric field 
component due to the gap voltage. 

The last boundary condition to enforce is that the 
magnetic field be continuous across the interface of Region 
I and Region II. Again using Equations 1-8, 1-10 and 1-12, 
and setting H $ (a,z) = H^ n (c,z), an expression for the phi 

component of the total magnetic field can be written as 





0 , for d < z < l 



for 0 <z <d 
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mm 



m =0 



W 2 




u n c n J l (u n a) cos 



nTtjz-h) 

. i J 



"= u {Eq. 1-16} 
Equation 1-16 can be reduced to one set of unknowns by solving 
for the coefficients b m analytically and using this result to 
solve for the c n in terms of b m and the unknown a n 's. This is 
accomplished by sifting out the desired coefficient using the 
orthogonal property of Fourier moment integrations which can 
be defined over the regions where the fields exist as 





(Eq. 1-17. a} 

0, m ^ n 

j, m = n #0 

q, m = n =0 



{Eq. 1-17. b) 
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{Eq. 1-17. c) 

Applying the integrations to both sides of Equation 1-15, the 
b m coefficients can be solved for in terms of known functions 
resulting in an expression for the terms in the second 
summation of Equation 1-16, as shown below: 




{Eq. 1-18} 



Therefore, 



( 2 ). 



V. bnH$‘Xv n a)l™.-y-£ll 



V 



O , (3) 



(Eq. 1-19} 



which can be rearranged to give 






V d j 0) m (2) / >. 

n l n \ v n a ) {Eq. 1-20} 



By applying the Fourier moment integrations to Equation 
1-16 then substituting in the expression for b m (Equation 
1-20) , the c n term can be written as 
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j 0 ( U n a ) 



N, 






(2) m 
in J 1 ( u n a ) k = 0 



X v k°kH} 2 \v k a)T kn 



where 



. v« v 1 'k 3> H ^ v k a K 

jae °-I L~— - T 



— 0 v k / t <0 // 0 <2, (v t a) 



ifc.n 



{Eq. 1-21} 



:.« = ^ COS ) cos ^-(z-/l)j^Z 




:n j 

~J 



r)-t) 



. [ kn h 
sin 



l 



(Eq. 1-22} 

We now have one equation with two unknowns. A second expres- 
sion for the unknown a n 's can be obtained by applying the 
moment integrations to the Equation 1-14 such that 



X 2 rr(2)/ . f lnnz\ Imz 

v n a n H o ( v „ a ) I COS^— jcos|— 

n = 0 J Q 

N 2 /• l 

= X U n c n J o(. u n<>) I i 

n = 0 J . 



dz 



( tin[z-h ) | //nTrzi , 
— — -I cos^—j—jdz 



{Eq. 1-23} 
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which gives 



V m a mH 0 



(2) 






(1) 



u n c n^o( u n a ) ^ m,n 

n= 0 



Rewriting Equation 1-21 yields 



(Eq.1-24) 






M-n^n'I 1 (ttn a ) I n ~ \c a k^l (Vk a )l'k,n 



it =0 



V 0 ^ 1 /t W H, ,2, (v k a) 



o V 

- JC0£ ot L 



(D u < 2 K 



To v * h' ' H 0 ( v k°) 



k,n 



{Eq. 1-25} 



Substituting this result into Equation 1-24 yields 



«. » SV. "a," " > 1 ~ l I 



« =0 J l(«n Q ) l Jt =0 






■/„ (u n a)T m , n J V 1 



; „?o" n / ®/ ,(«. - ) \ v. / t (1) H o (i, (v t . ) 



X- 1 



^ Jk.n 



{Eq. 1-26} 



Collecting terms and rewriting yields 
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{Eq. 1-27} 



This can be written in matrix form as 



N 



^ ^m,k a k ~ B / 



for m = 0, N j 



* =0 



{Eq. 1-28} 



where 



( 2 ) 



A m,k =v k H 1 (v k a)V mtk m *k 



A,», m =%"i (J V m a)P m , m - vJnf\v m a)J ™ 



V, V 1 t V,H ?’ ,v k a '> 

, p - 



w , (Dr, (2),„ * 

k = 0 V k^ k H 0 ( V k a ) 



{Eq. 1-29} 



and 






^ w /q(« n a) T m t n T k,n 



p = y 

m ,k ( 2 ) 

n=0 {£q. i-3 0 } 

The current as a function of distance along the antenna 
is then computed by 
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I (z ) = 2 k a H Q ( a,z ) 



for d <z <h 




n =0 



■)} 




{Eq. 1-31} 



The input impedance is then given as the ratio of the gap 
voltage Vo to the current at the end of the gap 



The computer program was written in FORTRAN 77 and was 
designed to run on an MS-DOS compatible system with a numeric 
coprocessor and 460 kilobytes of available memory after the 
disk operating system in loaded. If the numeric coprocessor 
is not available, the source code could be recompiled to 
emulate the specific coprocessor calls. The user would notice 
a dramatic increase in run time, but as all calculations are 
done in single precision, the accuracy should not be affected. 

The program code is divided into two separate programs. 
The first program, called MONO. FOR, is the main program that 




{Eg. 1-32} 



B. SOFTWARE OVERVIEW 
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calls a number of subroutines in SUBS. FOR. Variable names 
used in the program code closely follow those used in the 
theoretical development. Each variable is explicitly declared 
in each subroutine and array sizes are kept to a minimum to 
reduce required memory at run time. 

The first section of code in MONO. FOR is used to input 
the physical parameters using either an unformatted ASCII data 
file as shown in Appendix E or by entering the values from the 
keyboard. In each case, measurements may be entered in either 
wavelength or meters. If a data file is used, the first line 
must be either a "w" or "m", and may be either upper or lower 
case, to distinguish between wavelengths or meters. Unit 
number zero was used for all READ statements, which allows use 
of the program in a batch mode where all inputs and screen 
directed outputs can be redirected by standard DOS methods. 
However, the use of compilers other than Microsoft Fortran 
(Versions 4.01 or 4.1) may produce compiler errors when using 
this designation. The user should consult the specific 
compiler manual for details on unit designations. 

Before prompting the user to identify the output file 
name, the program checks for two conditions that will produce 
inaccurate results. The first condition is when the upper 
plate height is a multiple of one half wavelength. Physical- 
ly, this condition will allow resonant cavity modes between 
the plates, introducing instabilities in the numerical 
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solution. The mathematical basis for this can be easily 
shown. When the upper plate height is a multiple of half the 
wavelength, i/ n will become zero for certain values of n. As 
can be seen in Table I, the value of a Hankel function ap- 
proaches -jeo as the argument approaches zero. This test 
occurs in the main program and supplies earlier error detec- 
tion. The singularity is not checked for in the subroutine 
that computes the Hankel functions (HAN1) found in SUBS. FOR. 

The second condition that is tested is when the distance 
of the top ground plane above the antenna, designated q, is an 
integer multiple of the antenna height. This is similar to 
the previous case in that a resonant area is present between 
the antenna and the upper ground plane. However, in this 
case, the value of ti n , which appears only in the P mk matrix, 
becomes zero for multiples of the ratio h to q. Figure 3 
illustrates the behavior of the Bessel function for a range 
of arguments [Ref. 4]. The zero order Bessel function, J 0 , 
is finite for small arguments; however the first order Bessel 
function, J lf (found in the denominator of Equation 1-30) 
rapidly approaches zero as the argument approaches zero, 
causing the errors. 

In an effort to reduce computation time at the expense of 
memory requirements, several arrays are dynamically loaded so 
that the values may be easily extracted to fill the system 
matrix ^ k and the driving vector B m . These include the 
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Figure 3 Bessel Functions of the First Kind 

values for the Fourier moment integrations, the Bessel 
function values and the P mk matrix. The T matrix has dimen- 
sions of M x x N x and consists of terms from Equation 1-23. 
The I matrix has dimensions of M x x 3 where the column 
positions correspond to the three equations in Equation 1-18. 
The Ja and P mk matrices are loaded in the subroutine ALOAD and 
passed back to the main program. The Ja matrix is Mj x 4 in 
size where the column positions represent the following 
values: 
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Column Position 



Value 



1 



K * J 0 (X a) 



2 



J 1 (u n a) 



3 



v n* H X ( v n a ) 

(vj* H 0 (u n a) 



4 



For large values of N w N 2 and M x , about 60% of the total 
computation time is involved in filling the last temporary 
array, the P mk matrix. This array is dimensioned to N x x M x 
and uses values from all three arrays for each location. By 
filling these arrays first, the computation time was decreased 
by an average of 85%, however the memory required increased 
by 300%! 

Once the system matrix and driving vector are stored in 
memory, it is a simple task to solve for the unknown a n 's 
utilizing the back substitution method. First the matrix is 
upper triangulated by performing an L-U decomposition, with 
pivoting, on the matrix using the FACTOR subroutine. Once 
this is accomplished, the SOLVE subroutine calculates the 
eigenvalues by back substitution. 

To maintain the convention that an outwardly travelling 
wave has a negative imaginary exponential argument, the CSR 
function was included in SUBS. FOR. Since not all compilers 
would return a negative imaginary value from the square root 
algorithm when the argument is a negative real value, the CSR 
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function insures the —j component of the square root of a 
negative argument is returned. 

The subroutines that compute the Bessel and Hankel 
functions for complex arguments and the matrix solution 
algorithms were supplied by Professor Michael A. Morgan of the 
Naval Postgraduate School. The Bessel and Hankel subroutines 
utilize the direct power series method for small values (less 
than or equal to five for the Bessel subroutine and six for 
the Hankel subroutine) and Hankel' s asymptotic formula for 
larger values. 

The output is directed to a user specified file on any 
drive or subdirectory as long as the entire string is less 
than 25 characters. Longer strings will truncate the extra 
characters without warning. Printed output includes all input 
variables including the number of coefficients chosen for 
Regions I and II and for the gap voltage expansion. Length 
measurements are reported in both meters and fractions of a 
wavelength. The system eigenvalues are then printed followed 
by the position number, distance from the end of the gap and 
magnitude and phase of the computed currents, up to the upper 
plate height. Finally the input resistance and reactance is 
computed and printed to the file. Sample input and output 
files are listed in Appendix E. 
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III. VALIDATION 



In order to validate any computational algorithm, one 
must establish the accuracy obtained by comparing the computed 
results against either experimental data or another widely 
accepted software package that has been proven to be accurate. 
For this case, the logical choice was to use the Numerical 
Electromagnetics Code to calculate the input impedance and 
current distribution for a range of different physical 
structures and compare these results against those obtained 
by MONO. From this point, the program sensitivity will be 
examined for various ground plate spacings, antenna lengths 
and radii. Chapter 4 will investigate the effect that 
reducing the number of modal expansion terms has on this 
accuracy. 

A. THE NUMERICAL ELECTROMAGNETICS CODE 

The Numerical Electromagnetics Code (NEC) is widely 
accepted as an accurate method of computing the antenna 
currents, input impedance and radiation patterns of antennas 
that can be modeled using wire or patch structures. NEC 
calculates these values by solving integral equations for the 
currents imposed on the structure by a voltage source or an 
incident electromagnetic wave. The structure is best modeled 
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by straight wire segments that should be less than ten percent 
of the operating wavelength and conforms to the physical 
geometry using a best fit approach. [Ref. 6, p. 3] 

The user supplies NEC with a formatted data file that 
specifies the location of the structure in the rectangular 
coordinate system, the number of segments used to compute the 
results, the driving source location and output parameters and 
the existence and characteristics of a ground plane in the XY 
plane. Because of the infinite number of images produced by 
a pair of perfect ground planes, the user must also construct 
enough images with the proper driving voltages that will yield 
consistent results. NEC then calculates the segment length, 
and uses this value as the gap distance between the ground 
plane and the base of the antenna. The currents are calcu- 
lated at the midpoint of each segment vice at the ends as is 
done in MONO. 

NEC uses two approximations in the calculation of the 
electric field, the thin-wire kernel and the extended thin- 
wire kernel . The thin-wire kernel can be used when the 
segment length to radius ratio is greater than 8, while the 
extended wire kernel is accurate when this ratio is at least 
2 . Since NEC sets the gap distance equal to the segment 
length and accuracy is improved if the segment lengths 
(especially across a boundary) are of equal size, setting the 
gap distance equal to the antenna diameter forces the segment 
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length-to-radius ratio to be equal to 2. Additionally, this 
gap distance is consistent with the requirement for a constant 
electric field potential inside the gap region. Therefore, 
the extended wire kernel, which uses a uniform surface current 
along the segment length and assumes no variation of the 
currents along the <f> direction, will be used for all valida- 
tion runs. These approximations also support the basic 
assumption of the TM mode dominating the currents on the 
monopole. For the conditions listed above, NEC's accuracy is 
within 1% of experimental results. [Ref. 6] 

B. COMPARISONS MADE 

The first goal will be to determine the minimum number of 
reflections that must be assembled in the NEC data set to 
obtain consistent results for a quarter wave antenna struc- 
ture. From this data, an appropriate number of reflections 
will be used that represents consistent results while reducing 
the calculation time for all future runs. This will be 
followed by demonstrating consistent results for various 
ground plane spacings. The upper ground plane was introduced 
to establish a closed region. The upper plate should, if far 
enough away from the antenna, have little effect on the input 
impedance and antenna currents for distances that are not 
exact multiples of half the wavelength. Again a value of / 
will be selected that represents reliable results but reduces 
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the calculation time required for MONO. The antenna radius 
will then be adjusted over the range of 0.005 A to 0.02 A 
followed by variations in the antenna height for a selected 
radius. 

For all MONO calculations, the number of modes used in 
Regions I and II (N lf N 2 ) will always be set to 60, while the 
number of modal expansion terms used for the gap field (MJ 
will be set to two times the upper plate spacing to gap 
distance ratio. Chapter 4 will further investigate the 
accuracy of the results as a function of the number of modes 
used. A preliminary check for consistent results can be 
obtained by evaluating the differences in the input impedance 
results. But these values may differ slightly since NEC 
calculates the current at the midpoint of each segment instead 
of at the ends of the segment as is done in MONO. Addition- 
ally, since the first segment used in NEC is the gap distance, 
the input impedance is actually calculated for a point not on 
the monopole but in the center of the feed region. In 
comparison, MONO calculates the input impedance at the base 
of the antenna. The effect this has on the comparisons will 
be discussed later in this chapter. Additionally, these 
differences may be emphasized since the magnitude of the 
current at the end of the gap region can be quite small so 
that even small differences in the calculated values will 
yield relatively high differences. Therefore, to better 
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assess the accuracy of the results, the calculated current 
distribution along the antenna will be compared graphically 
with that obtained from NEC. 

C. RESULTS 

Since the ground planes used to divide the area into 
regions create an infinite number of reflections, one must 
first examine the minimum number of reflections required to 
obtain convergence to the half space case, having no upper 
ground plane. For this set of runs, the antenna is 0.24 A high 
with radius of 0.01A and the upper ground plate was placed 1.4 
wavelengths from the lower plate. Appendix F has a complete 
listing of some of the NEC data sets used for the validations. 
It should be noted that the distance between the lower plate 
and the center of the reflected dipole is equal to twice the 
plate spacing. This is due to the complete imaging, including 
the lower ground plane, of the physical structure from the 
upper ground plane. One can see from Figure 4 that adding a 
second reflection does not change the current distribution by 
any appreciable amount, therefore one reflection will be used 
for all further validations. 

Now that a baseline data set has been decided upon, the 
effects of various plate spacing, antenna heights and radius 
will be examined. The first item to look at will be the 
effect the plate spacing has on the current distribution and 
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Antenna Current Distribution (NEC) 

As the Number of Reflections Vary 
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Figure 4 Antenna Current vs Number of Reflections 
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Distance along Antenna (meters) 



the calculated input impedance. Table II shows the numerical 
results for an antenna of height 0.24A and radius of 0.01A 
placed between ground planes that are spaced 3.1, 1.4 and 0.81 

Table II RESULTS FOR VARIOUS UPPER PLATE SPACINGS 



Plate 


NEC 


MONO 


% Diff. 


Spacing 


Rin 


Xin 


Rin 


Xin 


R 


X 


3.10 A 


44.63 


13.35 


42.65 


13.14 


4.64 


1.60 


1.40 A 


44.90 


13.97 


43.14 


13.95 


4.08 


0.14 


0.81 A 


43.51 


14.84 


44.67 


13.84 


2.66 


6.74 



wavelengths. Figures 5 through 7 show the current distribu- 
tion for each of the above cases. While none of the results 
exceed acceptable limits, the case of /=1.4A is closer to the 
results computed by NEC and will be used as the baseline for 
further comparisons. It should be noted here that while the 
example values selected for the plate spacing have little 
effect on the result, the plate spacing can not be an integer 
multiple of a half wavelength. Likewise, the distance between 
the top of the antenna structure and the upper plate can not 
be a multiple of a half wavelength. Physically, these 
configurations would present the ability to support resonant 
modes that do not require a driving source. In reality, this 
condition will not exist as there is always a source of 
resistance which causes all modes to decay in time. However, 
the first condition will exist numerically where H 0 <2) ap- 
proaches negative infinity when »/ n (Equation 1-11) equals 
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NEC vs MONO 

Upper Plate Spacing = 3.1 Lamda 
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Figure 5 Comparison of NEC vs MONO when / =3 . 1A 
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Figure 6 Comparison of NEC vs MONO when / =1.4 A 
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Figure 7 



Comparison of NEC vs MONO when l =0.81A 



Distance along Antenna (meters) 



zero. By replacing the wave number k Q with 2n/\, then i/ n is 
zero when the plate spacing / is an integer multiple of a half 
wavelength. A similar condition will also occur when the 
distance between the top of the structure and the upper plate 
is a multiple of a half wavelength. In this case, J x found in 
the denominator of P m<k (Equation 1-30) is zero when is zero. 
The program code will alert the user to both conditions and 
halt execution. 

The next variable to consider is the antenna height. For 
extremely short antenna structures, the current decays almost 
linearly along the wire from its initial value to zero. When 
the antenna is about one quarter wavelength long, the current 
distribution takes on the familiar cosine shape. As the 
antenna length approaches one half wavelength, the distribu- 
tion approaches a sine wave shape and the input impedance 
approaches infinity as the current approaches zero. 

Table III lists the results when the structure is placed 
between the two ground planes spaced 1.4A apart for various 
antenna heights and radii. In the case where the antenna 
height approaches a half wavelength, the percent difference 
in the input impedance values appears to be large. This can 
be attributed to the differences in the location on the 
structure where the input impedance is actually calculated. 
This disparity combined with the low magnitudes in general can 
lead to alarming discrepancies with the NEC results. 
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Table III NEC, MONO RESULT COMPARISONS FOR VARIOUS h AND a 



Antenna ( x ) 


NEC 


Mono 


% Diff 


h 


a 


Rin 


Xin 


Rin 


Xin 


R 


X 


.12 


.005 


6.6 


-137.9 


7.4 


-148.7 


10.7 


7.3 


.24 


.005 


42.0 


12.4 


40.6 


10.8 


3.3 


14.8 


.48 


.005 


254.8 


-285.0 


339.9 


-281.7 


25.0 


1.2 


.12 


.01 


7.7 


-104.6 


9.3 


-118.2 


17.8 


11.5 


.24 


.01 


44.9 


14.0 


43.1 


14.0 


4.0 


0.1 


.48 


.01 


177.7 


-203.8 


264.8 


-199.9 


32.9 


2.0 


. 12 


.02 


9.8 


-76.3 


13.4 


-92.2 


26.5 


17.2 


.24 


.02 


47.9 


13.4 


46.1 


15.9 


4.0 


15.5 


.48 


.02 


129.6 


-139.8 


220.1 


-120.9 


41.1 


15.6 



Therefore, one roust look at the current distribution to 
accurately determine the differences in the computed results. 
Figures 8 through 11 show that the differences in the 
magnitude of the current distribution are within two percent. 
Additionally, one could also compare the phase difference 
between the two methods to further establish an acceptable 
method of judging the accuracy. Figure 12 shows that the 
difference between the phase of the current as calculated by 
MONO and NEC is well within two percent along the entire 
length of the antenna. 

Another point that can be observed is the effect that a 
change in the antenna radius has on the results. Once again 
a review of the current distribution is required to accurately 
assess the effect any errors would have on the far field 
pattern; however one detail is evident. As the antenna 
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becomes thinner the difference between the two results is 
reduced. This effect can be attributed to the increased 
resolution in the NEC results as the radii and, therefore, the 
gap distance are reduced. Additionally, the approximation to 
a TM current distribution along the surface of the structure 
becomes more accurate with a thinner structure. 

In review, it has been shown that due to the differences 
in the methods, the best method to assess the accuracy of MONO 
using NEC is obtained by comparing both the current magnitude 
and phase distribution along the antenna surface vice the 
input impedance values. Using this criterion, MONO is 
consistently within two percent of the results obtained with 
NEC for a variety of structures. Chapter 4 will investigate 
the difference caused by reducing the number of expansion 
terms and attempt to identify a minimum number of terms 
required to obtain results that are consistent with those 
obtained here. 



37 




(cojodujD) ©pn4|u6D^ juejjno 



Figure 8 NEC vs MONO for Quarter Wave Monopole 
(Radius = 0.005A) 
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Figure 9 NEC vs MONO for Quarter Wave Monopole 
(Radius = 0.02A) 
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Figure 10 NEC vs MONO for Antenna Height = 0.12 A 
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Figure 11 NEC vs MONO for Antenna Height = 0.48A 
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Figure 12 Phase Comparison between NEC and MONO 
( h = 0.48A, a = 0.02A) 
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IV. CONVERGENCE CONSIDERATIONS 



In any realizable computing machine where there is a 
finite quantity of storage, the first question that must be 
answered is the number of terms required to obtain accurate 
results. If one uses too few terms, the results are er- 
roneous. However, using too many terms can often lead to 
inaccuracies as well. If the additional terms are small 
enough to be approximately on the same order as the precision 
of the computer, then the additional terms may actually 
increase the error. To compound the problem, the additional 
terms require additional memory and increase the time required 
to compute the results. In investigating this question, we 
will first examine the number of terms required for Regions 
I and II (Nj and N 2 ) followed by an analysis of the number of 
terms required to represent the electric field in the gap 
region (M x ) . 

A. MODAL TRUNCATIONS FOR REGIONS I AND II 

In the previous chapter, the number of modes chosen for 
the two regions was the maximum allowed in the program, namely 
60. But is that enough? One could say that since the results 
compare favorably with NEC, the solution has converged. A 
more accurate approach is to increase the number of modes and 
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compare the results. However, due to the memory limit of 640 
kilobytes of addressable memory imposed by the compiler, the 
source code was recompiled on a 80386-specific FORTRAN 
compiler that allows for larger arrays and is only limited by 
the available memory of the computer being used. Figure 13 
compares the current distribution of a quarter wave antenna 
that has the upper plate height set at a large value (3.1A) 
for 60 modes and 100 modes. This additional information 
proves that the system has converged to an accurate solution 
with 60 modes and any additional modes are not required. The 
effect reducing this number has on the accuracy will be 
discussed later in this chapter. 

B. MODAL NEEDS FOR THE GAP REGION 

As mentioned earlier, the original approach was to expand 
the electric field along the surface of the antenna in one 
expansion between the ground planes as shown below 



E z (a,z) = 



{ 




0 , 



d ’ 



o 



for 0 <z <d 



for d <z <h 
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Modal Convergence for N1,N2 

(h=0.24. a = 0.01, d = 0.02, 1 = 3.1, M1=300) 
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Figure 13 Results for N x = N 2 = 60 and 100 Modes 



45 



Distance along Antenna (meters)* 



However, this approach failed to yield consistent results and 
a review of the source code showed no errors in the logical 
design. Upon an examination of the electric field in Region 
I, it was discovered that a very high value for N x would be 
required in order to obtain convergence to the known field in 
the narrow feed gap. It was believed that significant errors 
in the electric field would cause inaccuracies in the result. 
The errors in the results were thought to be caused by the 
lack of terms used to describe the electric field. The 
electric field in the gap region can be thought of as a step 
function of short duration. In the frequency domain, func- 
tions of this type occupy a large spectral bandwidth, requir- 
ing many terms for the description. The large expected value 
of N x would necessitate the inversion of a very large matrix, 
thus making this approach inefficient. To avoid this problem, 
a separate expansion for the constant electric field was added 
to Equation 4-1 to form Equation 1-14. 

Once again the number of terms required was the main 
question. However, this field is known prior to the computa- 
tion and the number of terms can be estimated using standard 
Fourier analysis. The Fourier transform of a step function 
is known to be a sine function with zero crossings at multi- 
ples of 1/d, where d is the duration of the pulse. Each lobe 
consists of impulse functions at intervals of the fundamental 
frequency, f ot which equals 1/T 0 , where T 0 is the period of the 
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pulse train. In this case, the upper ground plane causes a 
reflection of the original structure, but with opposite 
polarity, at a distance of 21 from the lower ground plane. 
Therefore the number of impulses in each lobe of the sine 
function equals 21/ d. This can be investigated graphically as 
well. Figures 14 to 19 (pages 48 to 50) illustrate the field 
computed for a structure with a gap distance ( d ) of 0.06 
meters and upper plate height of 9.54 meters and fed by a 
constant one volt source for various numbers of terms in the 
series. As can be expected, when the number of terms are not 
an integral multiple of l/d, excess oscillations occur. It is 
also important to note that as the number of terms increase, 
the calculated curve approaches the ideal value at the 
discontinuities, while the pointwise error along the rest of 
the curve increases. However, the error in the least square 
sense decreases. This is known as the Gibbs phenomenon. Near 
a discontinuity for a large number of terms, the overshoot is 
approximately 8.95 percent of the desired value, or -18.17 
volts [Ref. 5 p. 557]. This value is reached when the number 
of terms equals 3 l/d as seen in Figure 19, therefore addition- 
al terms will not reduce the error at the discontinuities. 
From these graphs, one can see that the number of terms used 
to represent the electric field in the gap region should be 
an integer multiple, less than or equal to three, of the ratio 
of upper plate height to gap distance. An upper plate height 
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E Field Convergence 

For 100 modes (0.631/d) 
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Figure 14 Gap Field Representation for 100 Terms 
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Figure 15 Gap Field Representation for 159 Terms 
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E Field Convergence 

For 300 modes (1.891/d) 
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Figure 16 Gap Field Representation for 3 00 Terms 
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Figure 17 Gap Field Representation for 318 Terms 
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E Field Convergence 

For 400 modes (2.521/d) 
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Figure 18 Gap Field Representation for 400 Terms 
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to gap distance ratio of 250 was considered to be the maximum 
desired. Therefore to insure the code would fit easily into 
the 640 kilobyte limit, the number of terms used for this 
expansion was set at 21/ d. This was the value used for all 
the validation calculations completed in Chapter 3. 

C. EFFECT OF REDUCING N, AND N 2 

If the number of terms can be reduced and still provide 
consistent results, one can then decrease the size of the 
internal arrays and reduce the computation time. For example, 
by reducing N x by 50%, the computation time is reduced by 
almost 48%, and when both N x and N 2 are reduced by 50%, the 
computation time is reduced by 71.5%! However, the current 
distribution along the antenna is unknown, making it difficult 
to calculate the number of coefficients required for consist- 
ent results based on the input parameters. The results 
obtained by independently varying N x and N 2 are shown in Table 
IV and are compared to the results obtained in Chapter 3 for 
an antenna with the following physical parameters: 



Antenna Height ( h ) 


0 . 24 A 


Antenna Radius (a) 


0.01A 


Gap Distance ( d ) 


0 . 02 A 


Plate Spacing (/) 


2 . 17 A 


Number of terms for 
Gap Region (Mi) 


217 
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Table IV 



COMPARISON OF RESULTS AS N x AND N 2 VARY 



N1 


N2 


Rin 


Xin 




Rin 


Xin 


60 


60 


43.08 


13.98 




N/A 


N/A 


50 


50 


42.70 


13.50 




0.88 


3.43 


40 


40 


42.80 


13.10 




0.65 


6.29 


30 


30 


41.60 


12.40 




3.44 


11.30 


20 


20 


41.20 


9.80 




4.36 


29.90 


50 


60 


42.20 


12.50 




2.04 


10.59 


40 


50 


41.80 


11.30 




2.97 


19.17 


30 


40 


39.90 


8.60 




7.38 


38.48 


20 


30 


35.60 


-0.80 




17.36 


105.72 


60 


50 


44.60 


16.70 




3.53 


19.46 


50 


40 


44.90 


17.47 




4.22 


24.96 


40 


30 


45.70 


18.99 




6.08 


35.84 


30 


20 


47.40 


22.60 




10.03 


61.66 



As expected, the difference in the input impedance increases 
as the number of terms is reduced. However, the interesting 
point is the effect that the apparent relationship between N x 
and N 2 has on the result. If N x and N 2 are equal, acceptable 
results are obtained for as low as 30 terms, while one must 
use at least 50 terms if the values are not equal, to get 
comparable results. A mathematical explanation for this 
relationship is not yet understood and will be investigated 
in future thesis efforts. Figure 20 represents the current 
distribution as N x and N 2 vary independently compared with the 
result when N x and N 2 equal 60. 
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Comparison of the Current Distribution 




Figure 2 0 Current Distribution as Nj and N 2 Vary 



D. EFFECT ON ACCURACY AS M, VARIES 

As seen in the previous section, one can reduce the 
memory requirements and computation time by reducing the 
number of coefficients used in Regions I and II. Using 
Fourier analysis, the user can calculate a value for the 
number of terms required to accurately represent the electric 
field along the gap (M x ) . Since this value is much larger 
than N x and N 2 , one could dramatically reduce the required 
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storage and speed up the computations by reducing Mj to the 
same order as N x and N 2 . 

To demonstrate the sensitivity of the input impedance to 
the number of terms used for the electric field along the gap, 
one can compare the results against those found when N : and N 2 
are 60 terms and is computed as 21/ d. Table V lists the 
results for same antenna structure described on page 51. As 
expected, the errors in the calculated input impedance begin 
to increase as the number of terms are reduced. However, this 
difference is not linear for all values of N x and N 2 . In fact. 
Figure 21 suggests an optimal number of M x terms, for given 
N : and N 2 , exists which results in smaller differences between 
the two cases. Since M x is a function of the gap distance, 



Table V MONO SENSITIVITY TO (M x ) FOR GIVEN Nj AND N 2 















k Dli 


N1 


N2 


Ml 


R 


X 




R 


X 


60 


60 


217 


43.08 


13.98 




N/A 


N/A 


60 


60 


60 


43.31 


13.87 




0.53 


0.77 


60 


60 


50 


43.02 


13.85 




0.14 


0.92 


60 


60 


40 


43.39 


13.95 




0.72 


0.20 


60 


60 


30 


42.44 


14.14 




1.49 


1.16 


50 


50 


60 


42.81 


13.44 




0.63 


3.85 


50 


50 


50 


42.59 


13.42 




1.14 


3.99 


50 


50 


40 


43.12 


13.53 




0.09 


3.21 


50 


50 


30 


42.13 


13.72 




2.21 


1.85 


30 


30 


40 


41.72 


12.39 




3 . 16 


11.36 


30 


30 


30 


40.70 


12.64 




5.52 


9.57 


30 


30 


20 


42.02 


13.26 




2.46 


5.14 


30 


30 


15 


39.41 


13.61 




8.52 


2.63 


30 


30 


10 


38.29 


14.73 




11.12 


5.38 
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the input impedances were recalculated when the gap distance 
is reduced to 0.01A for selected values of N x and N 2 as shown 
in Figure 22. While this clearly demonstrates that an optimal 
number of terms may exist that produces results which approach 
those for large values of M lf N x and N 2 , a clear relationship 
between the physical structure and these values are not yet 
fully understood. 

But how can the number of terms be reduced by almost 91 
percent, and yet the difference is no greater than 5.13 
percent? First, the basis for selecting the number of terms 
used was based on obtaining accurate representation for the 
electric field along the antenna. However, since the far 
field patterns are derived directly from the currents, it is 
the current that must remain consistent and not necessarily 
the electric field along the surface of the antenna. From 
Equation 1-31, we see that the current is derived from the H^, 
or the first derivative of the potential field with respect 
to p. However, the electric field along the z direction is a 
function of the second derivative of the potential with 
respect to z. Since the derivative process increases the 
noise, or errors, it is evident that small differences in the 
potential field may yield wide variations in the electric 
field while the magnetic field, and therefore the current 
distribution, may remain within accepted accuracy limits. To 
prove this point, Equation 4-2 shows an expression that 
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Average Error vs Ml 

For given number of terms N1 <= N2) 




N1=60 
N 1=50 
N 1 =40 
N1=30 



Number of Modes Used for Gap Region (Ml) 

Figure 21 Average Differences for Reduced Values of M x for 
Given N x and N 2 Values. (Gap = 0.02 A) 



Average Error vs Ml 

For given number of terms N1 (= N2) 



N 1 =60 
N1=40 
N 1 =30 



10 20 30 40 50 60 

Number of Modes Used for Gap Region (Ml) 

Figure 22 Average Differences for Reduced Values of M x for 
Given N x and N 2 . (Gap = 0.01A) 
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approximates the electric field along the surface of the 
antenna by taking the second derivative of Equation 1-31 with 
respect to z. 




+ jCO£ 



V 



A/, 



m = 0 



,3, <M 



m 



m * m 



COS 



\t1lTlZ 

T~, 



{Eq. 4-2} 

This equation can now be used to visualize the approximate 
change in the electric field along the antenna for a given 
number of modes Mi. As shown in Figure 23, the change in the 
magnitude of the current distribution when is changed from 
217 (which equals 2 I/d) to 30 is negligible. However, Figure 
24 shows the curves obtained from Equation 4-2 for the same 
parameters as used in Figure 23. As one can see, the field 
is accurately represented near the end of the wire, however 
near the gap region, there are an insufficient number of 
coefficients to accurately compute the fields in this region. 
This refutes the original assumption that the electric fields 
must be accurately represented to guarantee an accurate 
current solution. Additionally, since can be on the same 
order as the values used for N x and N 2 , then the coefficients 
used to describe the gap field and the field along the 
structure can be combined into one coefficient. The next 
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chapter will explore the theory and the differences this new 



approach has 
coefficients 



over calculating a separate set of expansion 
for the gap region. 
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Current Distribution For a Given Ml 

For N 1 =N2=60 




Distance along Antenna (meters) 

Ml =30 Terms Ml =2 17 Terms 



Figure 23 Current Distribution as Ml Varies (h =0.24 A) 
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Figure 24 Approximated Electric Field Distribution 
(Calculated using Equation 4-2) 




59 



V. A SIMPLER APPROACH 



During the early stages of the software development, it 
was decided that a separate expansion for the constant 
electric field in the gap region was required to obtain 
consistent results for the input impedance and current 
distribution along the structure. This decision was based on 
the belief that large errors in the total electric field along 
the structure would result in large errors in the current 
distribution. However since the current distribution is 
derived directly from the magnetic field, and is relatively 
insensitive to the electric field, the number of terms used 
for the gap region expansion was reduced to the same order as 
those used for N* and N 2 without introducing appreciable 
errors. Noting this, the original method of describing the 
system was reevaluated such that one could reduce the required 
memory and speed up the calculations even further by combining 
the flj, and b n coefficients used in Equation 1-10 (page 11) into 
one set of terms and by deleting the separate expansion along 
the gap region entirely. This chapter will work backwards 
from the technique presented in Chapter 2 and will note the 
changes in the system equations resulting from this approach. 
The resulting code, which will be referred to as M0N08, will 
then be compared to the results obtained for the cases with 
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the largest differences between NEC and the code using the 
separate expansion (called M0N07) . 



A. MODIFYING THE TECHNIQUE 

As was discussed in Chapter 2, the previous approach used 
two expansions in Region I where p > a. The first expansion 
sets the tangential electric field to zero from the lower 
ground plane to the top of the antenna. A second expansion 
sets the tangential electric field to - V 0 /d between the lower 
ground plane and the feed point of the antenna and sets E z to 
zero from this point to the upper ground plane. Since the 
Fourier transform of a sum is the sum of the transforms [Ref. 
7], we can combine the expansions for the potential field in 
Region I as shown below 



where a n equals the sum of the a n and b m coefficients used in 
Chapter 2. Again we match the magnetic fields across the 
boundary between Region I and II such that 




{Eq. 5-1} 




n =o 
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By utilizing the orthogonality property of the Fourier moment 
integrations used in Equation 1-15, the system equation can 
be expressed as 



N i 

I 

1 = 0 



N 



2) ( v * a )\ X 



u n^0 ( ll n a ^ T m ,nT k,n ^ 2 uP^(\> n\ 1 ^ 

^ / ‘ V m a m tl 0 



„=0 Cl 



V _0, ( 3 ) 
d 



~ y ~7T >n 



{Eq. 5-3} 

for m = 0 to Nj. This can be expressed in matrix form as 
N i 

y, A m,k a k = for m =0,/^! 

1=0 

(Eq. 5-4} 

where A„, k is given in Equation 1-29. As the b m coefficients 
are now included in the a n terms, the driving vector, B m , is 
reduced to 



= y I m 



(3 )V e 



{Eq. 5-5} 

Once the unknown a n coefficients are found, the current along 
the antenna, as a function of z, is now given by 



j 

1 (z) = 2® ^ < 2> (v„ a )cos(^) 

n = 0 

{Eq. 5-6} 

B. SOFTWARE MODIFICATIONS 

As seen in Equations 5-3 to 5-6 above, removing the 
separate expansion for the electric field along the gap region 
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requires changing the software code for only the driving 
vector B m and the expression for the current as a function of 
position along the antenna structure. The first advantage 
realized was the reduction in the amount of memory storage 
required to run the program. Since it was shown in the 
previous chapter that the system had converged when the number 
of modes in Regions I and II were set at 60, independent of 
the number of modes in the gap region, it was not necessary 
to dimension any array greater than 60. However, we can now 
reduce the size of the I, T, Ja, and P. k matrices by almost an 
order of magnitude, such that the amount of memory required 
to run the program has decreased by 75 percent to less than 
104 kilobytes. Additionally, the execution time of the 
program has decreased by an average of 35 percent for the runs 
presented here. 

But does this method yield results that are comparable to 
the results when the separate expansion is used? We will 
examine the differences between M0N07 and M0N08 for various 
plate spacings, antenna heights, as well as the effect that 
independently varying N x and N 2 has on the results obtained 
earlier. 

C. RESULTS 

As a first test of the accuracy for M0N08, we will 
examine the differences in the calculated input impedance 
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between M0N07 and M0N08 as the ground plane spacing is varied. 
The three structures used in Chapter 3 for Table II were rerun 
where N x and N 2 were set to 60. Table VI indicates that the 

Table VI COMPARISON OF M0N07, M0N08 FOR VARIOUS / 



Plate 


M0N07 


M0N08 


% Diff 


Spacing 


Rin 


Xin 


Rin 


Xin 


R 


X 


3.10A 


42.65 


13.14 


43.09 


13.12 


1.03 


0.15 


1.40A 


43.14 


13.95 


43.15 


13.96 


0.02 


0.07 


0.81A 


44.67 


13.84 


44.68 


13.82 


0.02 


0.14 



results obtained from the new method differ only slightly from 
the results obtained when the gap expansion is used. 

The next test of the accuracy considered different 
antenna heights. In this case, the comparison is made between 
the structures that resulted in the highest difference between 
M0N07 and NEC, or when the antenna radius was 0.02 A. Table 
VII shows the results of this comparison for a plate spacing 
of 1 . 4 A . Once again, the calculated input impedances differ 
very little from those obtained using the separate gap 
expansion. 



Table VII COMPARISON OF MONO 7 , M0N08 FOR SELECTED h 



Antenna 


M0N07 


M0N08 


% Diff 


Length 


Rin 


Xin 


Rin 


Xin 


R 


X 


0 . 12 A 


13.35 


-92.19 


13.32 


-92.13 


0.22 


0.06 


0 . 24 A 


46.08 


15.90 


46.08 


15.87 


0.00 


0.18 


0 . 48 A 


220.12 


-120.91 


219.44 


-121.34 


0.31 


0.35 
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As noted in Chapter 4, the results are more consistent 
with the solutions obtained from NEC when N x equals N 2 . Figure 
25 illustrates that this same phenomenon occurs for M0N08 . 
Further investigation and study is required to understand the 
relationship that the number of modes in Regions I and II have 
on each other. Once this is understood, it may be possible 



Comparison of the Current Distribution 




Figure 25 Current Distribution using M0N08 as N x and N 2 Vary 

to derive an expression for the minimum number of modes 
required to obtain results within a specified accuracy range. 
This expression could then be used to allow the software to 
dynamically adjust the number of modes required to obtain 
accurate results at run time. 
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In all the above cases, the results were compared to the 
solution using M0N07 where Mi was chosen as 21/ d. However, 
Figures 21 and 22 (page 56) indicates that there is an optimum 
number of modes for the gap region to reduce the differences 
between the solutions of MONO and NEC. Table VIII shows that 
M0N08 returns the same result for the input impedance as M0N07 
did where M x , N x and N 2 were equal. While this result may not 
be the closest to NEC's results, the differences are consis- 

Table VIII COMPARISON OF M0N08 AND M0N07 AS M x VARIES 
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13.42 


13.44 


13.42 


13.53 


13.72 


30 


40.70 


41.61 


41.66 


41.72 


40.70 




12.64 


12.32 


12.32 


12.39 


12.64 



tently less than ten percent. 

In conclusion, consistent results are obtained without 
the separate expansion for the gap region. While these 
results may not offer the closest possible solution to that 
available by NEC, the small difference in the solution 
obtained without the extra terms is more than offset by the 
significant reduction in required memory and calculation time. 



66 



When the relationship between the physical structure and the 
minimum number of modes that are required to obtain the best 
possible solution is understood, then the software can be 
modified so as to select the number of terms required at run 
time. Until then, M0N08 should be the method of choice as 
its results are more predictable and is capable of returning 
the result within acceptable accuracy limits without "trial 
and error" . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

We have shown that one can calculate the current dis- 
tribution along a monopole using multiregional cylindrical 
harmonic expansions that are consistent with the results 
obtained from the Numerical Electromagnetics Code (NEC) . 
While this software package is not intended as a replacement 
for NEC, it is the beginning step in the development of 
software directed towards solid antenna structures that are 
rotated about a symmetrical axis with a simple and flexible 
input structure. The particular software that has been 
developed has the following characteristics: 

• Simple input data structure. All inputs entered direct- 
ly, via unformatted ASCII file, or redirected using 
standard DOS redirection codes allowing batch processing 
of multiple runs. 

• Very flexible for a given monopole antenna structure. 

• Can be readily modified to allow for a homogeneous 
dielectric coating on the surface of the antenna. This 
coating could be used to modify the surface currents and 
the far field radiation pattern. 



As currently written for the simple monopole, the user is 
able to modify all the physical parameters of the structure 
including driving source frequency. With some modifications, 
including the incorporation of a finite element algorithm in 
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one sub-region, the user may also include the effect of a 
nonhomogeneous dielectric coating so that the far field 
radiation pattern may be tailored as desired. 

One difference that should be noted between the techni- 
ques used by NEC and MONO is the location used to calculate 
the input impedance. Due to the NEC's segmentation of the 
entire structure, the input impedance is actually calculated 
in the center of the gap region. However, MONO calculates the 
input impedance based on the value of the current at the base 
of the antenna. This difference may lead to large discrepan- 
cies in the calculated input impedances if the value of the 
current is changing rapidly at the end of the gap region, as 
was observed for monopoles whose height approaches a half 
wavelength. Other differences include: 

• MONO requires longer computation time for larger struc- 
tures than does NEC. 

• MONO Requires larger available memory. 

• MONO currently has limited error detection and correction 
within the subroutines. 

• The current inability to calculate a number of internal 
parameters requires a "trial and error" approach to 
determining the minimum number of modes required for 
proper convergence. 

• User is required to make extensive modifications to 
MONO's code for various types of antenna structures. 

• MONO is a highly specific, concept validation platform 
vice a generally applicable software package. 
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B. RECOMMENDATIONS 

Continued research is required in order to understand the 
effect that the apparent relationship between N x and N 2 has on 
the accuracy of the result. Additionally, there appears to 
be an optimal number of terms required for the gap region that 
results in an apparent improvement of the result when compared 
to those obtained using NEC. These unanswered questions point 
out the need for additional analysis of the relationship 
between the number of terms required to represent the electric 
field in a specific region and the overall accuracy of the 
final result. When this is accomplished, modifications should 
be made to allow the software to calculate the optimal number 
of terms required for each expansion that would yield results 
in the shortest time possible for a specified accuracy range. 

Another modification could be made to the code that would 
reduce the amount of memory required to run the program. This 
can be accomplished by temporarily storing only one row of the 
P m k matrix at a time, vice storing the entire matrix. This 
approach could reduce the amount of required memory by as much 
as 50 percent. 

To better assess the benefits of this method over the one 
used in other computational schemes, experimental data should 
be compared to the computed results for a variety of struc- 
tures. In this way, a basic understanding of when this method 
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best approximates the actual performance characteristics of 
an operational antenna can be obtained. Once these areas are 
better understood, a collection of routines for various 
physical structures should be developed using the basic 
technique developed in Chapter 2. This collection could 
include top loaded monopoles, and conical and other structures 
that feature symmetry about a axis of rotation. 
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APPENDIX A 



SOFTWARE FLOW CHART FOR MONO WITH GAP FIELD EXPANSION 
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ALOAD 

Subroutine 
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APPENDIX B 



SOURCE CODE FOR MONO WITH GAP FIELD EXPANSION 



PROGRAM MONO 
C (Revision 7.3) 

C 

C Computing Currents on a Monopole Between Parallel 
C Plates using Cylindrical Harmonic Expansions in 2 Regions. 

C Links with SUBS4.0BJ. Created 05/87 LT. R.C. Hurley 

C 

Variables 

= Frequency in MHz 
= Monopole height in Meters 
= Monopole radius in Meters 

= Distance between lower ground plane and feed point 
in Meters 

= Distance between ground planes in Meters 
= height from top of monopole to upper ground plane 
in Meters 
= Imaginary One 
= Value of pi (3.1415927) 

= 2*pi*F/c 

F = frequency of operation in Hz 
c = Speed of Light 
= K 0**2 

= Number of Coefficients in each of the two regions 
= Size of the arrays used in FACTOR and SOLVE 

= Screen Output 

= Main Output File to Disk (Name selected by User.) 
***** Variable Definitions ***** 

INTEGER k,m,M1,Mmax,mpn,n,ns,Nl,N2,Nbig,NI,Nmax,P(61) 

REAL a,ai,Am,An,Bl,c,CI ,CM,CP,CR,d,di ,DeltaZ,Denom,eps, FMHZ,h,hi 
REAL I (0:500, 1 :3),K0,K2, l, li , Lambda, Pi ,q,Rin,S 
REAL T(0:500,0:60),w,Xin,y,yhat 

COMPLEX Amk(61 ,61 ),B(61 ),C1 ,C2, IZ(61 ), j, JO, J1 , Ja(0:500, 1:4) 
COMPLEX Pmk( 0:60, 0:500), Sum, Suml ,Y0,Y1 ,Zin 
CHARACTER*25 Fname 
CHARACTERS Bell 
CHARACTERS HOW 
CHARACTERS INPUT 
CHARACTERS2 INname 
COMMON/Blk1/I,j,Nl,Pi,d,Nbig 
COMMON/Blk2/K2,l,Mmax,N2,q,a,T 



FMHZ 

h 

a 

d 

l 

q 

j 

pi 

K0 



where 



K2 

N1,N2 

NBIG 



Unit 

Unit 



********************************* 

* Main Program Code * 

********************************* 



C 

C 

C 

Bel l=CHAR(7) 



*********************** 

* Input Section * 
*********************** 
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Nbig=61 

c=3.0E+08 

WRITE(6,*) 'Do you wish to use an input data file [Y/N]?' 

READ(0, 100) HOW 

IF (HOW.EQ. 'y # .OR.HOU.EQ. ' Y # ) THEN 

WRITEC6,*) 'Enter input data file name with extension 
READ(0, 100) INname 
OPEN(Unit=4, Fi le= INname) 

READ(4, 100) INPUT 

IF ((INPUT .EQ. 'w') .OR. (INPUT .EQ. 'W') .OR. (INPUT .EQ. 

+ 'm') .OR. (INPUT .EQ. 'M')) THEN 

READ(4,*) FMHZ,hi ,ai ,di , l i ,N1 ,N2,M1 ,NI 

ELSE 

WRITE(*,*) Bell 

WRITE(6,*)'Error detected in first line of data file ' 

+, INname 

WRITE(6,*) 'First line must be either M or W for meter 
+ values or wavelength values.' 

STOP 

ENDIF 

ELSE 

2 WRITE(6,*) 'Enter data in either [M]eters or [Wlavelength [H 

+/W]?' 

READ(0,100) INPUT 

IF (INPUT. EQ.'m'. OR. INPUT. EQ.'H') THEN 

WRITE(6,*) 'Enter Frequency in MHz: ' 

READ(0,*) FMHZ 

WRITE(6,*) 'Enter Monopole Height, h, in meters: ' 

READ(0,*) hi 

WRITE(6,*) 'Enter Monopole Radius, in meters: ' 

READ(0,*) ai 

WRITE(6,*) 'Enter Feed Gap, d, in meters: ' 

READ(0,*) di 

WRITE(6,*) 'Enter Upper Plane Height, l (meters): ' 

READ(0,*) li 

ELSEIF (INPUT. EQ.'W'. OR. INPUT. EQ.'W') THEN 

WRITE(6,*) 'Enter Frequency in MHz: ' 

READ(0,*) FMHZ 

WRITE(6, # ) 'Enter Monopole Height, h, as factor of wav 

+elength : ' 

READ(0,*) hi 

WRITE(6,*) 'Enter Monopole Radius, as factor of wavele 

+ngth : ' 

READ(0,*) ai 

WRITE(6,*) 'Enter Feed Gap, d, as factor of wavelength 

+ : ' 

READ(0,*) di 

WRITE(6,*) 'Enter Upper Plane Height, l, as factor of 
♦wavelength : ' 

READ(0,*) li 

ELSE 

WRITE(*,*) Bell 

WRITE(6,*) 'You must enter either an M or W. 

♦Reenter.. ' 

GOTO 2 

ENDIF 

WRI TE(6,*) 'Enter # Region 1 Coeffs. (N1 .LE. 60): ' 

READ(0,*) N1 

WRITE(6,*) 'Enter # Region 2 Coeffs. (N2 .LE. 60): ' 

READ(0,*) N2 

WRI TE(6,*) 'Enter # Gap Region Coeffs. (Ml .LE. 500): ' 

READ(0,*) Ml 

WRITE(6,*) 'Enter No. Points for I(z) including ends: ' 

READ(0,*) HI 
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ENDIF 

IF (M1.GT.500) Then 
WRITE(\*> Bell 

WRITE(6,*) 'Ratio l/d too large. Adjust either:' 
WRITE(6,*) ' a) l or d such that 2l/d < 500' 

WRITE(6,*) ' b) Array sizes of Im,T,Ja' 

STOP 

ENDIF 

*********************************************** 
* * 

* Calculate the Physical Parameters * 

* in both Uavelenght and Meters * 

* * 
*********************************************** 



Lambda=c/(FMHZ*10**(6)) 

IF (INPUT. EQ.'w'. OR. INPUT. EQ.'W') THEN 
h=hi*Lambda 
a=ai* Lambda 
d=di*Lambda 
1=1 i*Lambda 

ELSE 

h=hi 

hi=h/Lambda 

a=ai 

ai=a/Lambda 

d=di 

di =d/ Lambda 
1 = 1 1 

l i= l/Lambda 

ENDIF 

*********************************************** 
* * 

* Is Upper Plate Height an * 

* Even Multiple of Wavelength? * 

* * 

*********************************************** 



IF (M00((l*2), Lambda). LT. IE-05) THEN 
WRITE(*,*> Bell 
URITEC6,*) ' ' 

WRITEC6,*) 'ERROR! Upper Plate Height is almost an even mul 
+tiple of half the wavelength.' 

WRITE (6,*) 'Ensure l is NOT within 0.00001 of an integer mul 
♦tiple of half the wavelength.' 

STOP 

ENDIF 

*********************************************** 

* * 

* Is Upper Plate Height an * 

* Even Multiple of Antenna Height? * 

C * * 

C *********************************************** 

c E=(10*h)- (ANINT((10*h)/l )*l ) 

IF ((M00(l,h)).LT. IE-05) THEN 
WRITE(*,*) Bell 
WRITE(6,*) ' ' 

WRITE(6,*) 'ERROR! Upper Plate Height is almost an even mul 
♦tiple of the antenna height.' 

WRITE(6,*) 'Ensure l is NOT within 0.00001 of an integer mul 
+tiple of h (in wavelength).' 



76 



STOP 



ENDIF 

WRITE(6,*) 'Enter output file name and extension: ' 

READ(0, 100) Fname 

OPEN (Uni t=7, Fi le=Fname) 

URITE<7,210) FHHZ 
WRITE(7,220) h,hi 
WRITE<7,230) a,ai 
WRITE(7,240) d # di 
WRITE(7,250) l,li 

WRITE(7,*) 'Region I Coeffs. (Nl): ',N1 
WRITE(7,*) 'Region II Coeffs. <N2): ' # N2 
WRITE(7,*) 'No. Gap Coeffs. (Ml): ',M1 
WRITE(7,*) 'No. Points for I(z): \NI 

j=(0. 0,1.0) 

Pi=3. 1415927 
K0=Pi*FMHZ/150.0 
K2=K0*K0 
q=l-h 

Nmax=MAX0(Nl,N2) 

Mmax=MAX0(Nmax,M1) 

URITE(6,*) 'Loading the Moment Matrices ' 

WRITE(6,*) ' ' 



C 

C 

C 

C 

C 

C 



*********************************************** 
* * 

* Evaluating Moment Integrations * 

* Over Region h < z < l. * 

* * 

*********************************************** 
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C 

C 

C 

C 

C 

C 



T(0,0)=q 
DO 11 k=1 ,N2 

T(0,k)=0.0 
CONTINUE 
DO 22 m=1,Mmax 
Am=m*Pi/l 

DO 22 n=0,N2 

An=n*Pi/q 

I F( (ABS(Am-An) ) .GT. 1 .OE-3) THEN 
D enom=Am* Am- An* A n 
S=SIN(Am*h) 

T ( m , n ) = - Am* S/ D enom 

ELSE 



mpn=m+n 

T(m,n)=0.5*q*((-1)**mpn) 

ENDIF 

CONTINUE 

URITE(6,*) ' Moment Matrix T Loaded.' 

WRITE(6,*) ' ' 



*********************************************** 



* * 

* Loading the Im Matrix for the * 

* Largest Number of Coefficients * 

* * 



*********************************************** 



DO 33 m=0,Mmax 

I F(m.EQ.O) Then 
I (m, 1 )=l 
I(m,2)=q 
I(m,3)=d 

ELSE 
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l(m,1)=l/2.0 

I(m,2)=q/2.0 

Am=m*Pi/l 

I(m,3)=SIN(Am*d)/Am 

END IF 

33 CONTINUE 

WRITE(6,*) ' Moment Matrix I Loaded.' 

WRITE(6,*) ' ' 

WRITE(6,*) 'Loading System Matrix ' 

CALL ALOAD(Ja,Pmk,Amk) 

URITE(6,*) 'System Matrix Amk Loaded.' 

WRITE(6,*) ' ' 

WRITE(6,*) 'Loading Driving Vector ' 

CALL BLOAD(FMHZ,M1,Pmk,Ja,B) 

WRITE(6,*) 'Driving Vector B Loaded.' 

WRITE(6,*) ' ' 

WRITE(6,*) 'Solving the Linear System .. Please Wait .... ' 
Ns=N1+1 

CALL F ACTOR ( Amk, P, Ns, Nbig) 

CALL SOLVE(Amk,B,P,Ns,Nbig) 

WRITE(7,*) 'Listing the Expansion Coefficients' 

URITE(6,*) ' Listing the Expansion Coefficients ' 

WRITE(7,*) 'N and B(N): ' 

DO 44 N=0,N1 

WRITE(7,*) N,B(N+1 ) 

44 CONTINUE 

Q *********************************************** 

C * * 

C * Computing Current from the * 

C * Expansion Coefficients * 

C * * 

Q *•****•**••••***•••••••••****•••••••••••••••••• 

WRITE(6,*) ' ' 

WRITE(6,*) 'Computing current on the monopole ' 

WRITE(7,*) 'Listing the Current on the Monopole' 
DeltaZ=(h-d)/(NI-1) 

WRITE(7,*) 'Position Distance I ( z) (A)' 

WRITE(7,*) ' Nuifcer (Meters) Mag Phase' 

eps=8.854E-12 

w=2.0E6*Pi*FMHZ 

yhat=w*eps 

y=d 

k=1 

55 Sum=(0.,0.) 

DO 66 n=0,Nl 

SuiPSum+ ( B( n+ 1 ) * J a ( n , 3 )*cos ( n*P i *y/ 1 ) ) 

66 CONTINUE 

Sum1=(0. 0,0.0) 

DO 77 m=0,M1 

Sun1=sum1+((I(m,3)*Ja(m,3)*COS(m*pi*y/l))/( Ja(m,4)*I(m,1))) 
77 CONTINUE 

IZ(k)=2.0*Pi*a*(Sum-j*yhat*Sum1/d) 

CM=CABS(IZ(k)) 

CR=REAL(IZ(k)) 

CI=AIMAG(IZ(k)) 

CP=180.*ATAN2(CI ,CR)/pi 
WRITE(7,200) k,y,CM,CP 
y=y+DeltaZ 
k=k+1 

IF(y.LE.(h+1.E-6)) GO TO 55 
WRITE(6,*) 'Complete.' 
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*********************************************** 
* * 



c 
c 

C * Computing Input Impedance * 

C * * 

C *********************************************** 

Zin=(1.0,0.0)/IZ(1) 

Rin=REAL(Zin) 

Xin=AIKAG(Zin) 

WRITE(7,*) 'Input Resistance = ',Rin 
WRITE(7,*) 'Input Reactance = ',Xin 

r *********************************************** 

C * * 

C * I/O Format Statements * 

c * * 

c *********************************************** 



100 FORMAT (A) 

102 FORMAK 1 5) 

103 FORMAKE12.3) 

200 FORMAK2X, I3,5X, F8.3,5X, 1PE12.3,5X,0PF8.2, 1PE12.3,5X,0PF8.2, 
+1PE12.3,5X,0PF8.2) 

210 FORMAT(' Driving Frequency, ',9X,F7.3,' MHz') 

220 FORMAK' Monopole Height, ',9X,F7.3,' meters (',F5.3,' Lambda)') 

230 FORMAK' Monopole Radius, ',9X,F7.3,' meters (',F5.3,' Lambda)') 

240 FORMAT(' Feed Gap Distance, ',9X,F7.3,' meters (',F5.3,' Lambda)') 

250 FORMAK' Upper Plane Height, ' ,9X,F7. 3, ' meters <',F5.3,' Lambda)') 

STOP 
END 
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Subroutine Source Code for MONO with Gap Field Expansion 



SUBROUTINE A LOAD ( J a, Pink, Amk) 

C (Revision 7.3) 

******************************************** 
* * 

* Variable Definitions * 

* * 

******************************************** 



C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



l = Distance between ground planes 

Radius = Radius of monopole 

q = Distance between top of monopole and upper plane 
d = Antenna feed gap 

K2 = K0**2 (From main program) 

T(m,n) = Moment integration result as function of m and n 

I(m,x) = Moment integration result of order x given m 

j = Imaginary One 

z = Intermediate value for Bessel Subroutine 

un, vn = Eigenvalues 

JO = Bessel Function of the First Kind of Order 0 

J1 = Bessel Function of the First Kind of Order 1 

HO = Hankel Function of the Second Kind of Order 0 

HI = Hankel Function of the Second Kind of Order 1 

SUM = Summation over n for the Amk matrix 

SI = Individual term used in SUM 

m = Column position of matrix Amk 

k = Row position of matrix Amk 

Hi = Number of subdivisons in Region 1 

N2 = Number of subdivisons in Region 2 

Amk = Complex Matrix value that describes the system 

This SUBROUTINE loads the complex A matrix that describes the system 
of size Nl by Nl. Maximum size of matrix is set at 50 by 50. This 
SUBROUTINE calles to the following subroutines found in SUBS. FOR: 

CSR - Complex square root that preserves negative imaginary 
portion. 

BES1 - Computes the Bessel Function values for a given z. 

HAN1 - Computes the Hankel Function values for a given z. 

Written 5/88 by R.C. Hurley using Microsoft Fortran 4.01. 

INTEGER k,m,Mmax,n,Nbig,Nl ,N2 

REAL l,a,q,d,K2,T(0:500,0:60),I(0:500,1:3),Pi,Rl,R2 
COMPLEX j,zv,zu,vn,un,J0n,J1n,H0n,H1n,SUM,S1 
COMPLEX CSR,Amk(Nbig,Nbig), Ja( 0:500, 1 :4) 

COMPLEX Pmk(0:60, 0:500) 
character’ll cr 

COMMON/Blkl/I , j ,N1 , Pi ,d,Nbig 
COMMON/Blk2/K2, l ,Mmax,N2,q,a,T 
cr=char(13) 

This section loads a matrix consisting of Bessel fuctions and 
other values required to load the Amk matrix. The colunns are 
filled as follows: 

Ja(n,1) = un*Jo(un*a) in ncmerator of Sum 

Ja(n,2) = J1(un*a) in denominator of Sum 

Ja(n,3) = vn*H1(vn*a) used for all locations 

Ja(n,4) = vn*vn*Ho(vn*a) used in diagonal elements 

WRITE(6 f 100) Mmax 
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DO 11 n=0,Hmax 

R1=K2-((n*Pi/q)*(n*Pi/q)) 

R2=K2-((n*Pi/l)*(n*Pi/l)) 

un=CSR(R1 ) 

vn=CSR<R2) 

zu=un*a 

zv=vn*a 

Call BES1(zu, JOn, Jin) 

Call HAN1 (zv, H0n,H1n) 

Ja(n,1) = Un*J0n 
Ja(n,2) = Jin 
Ja(n,3) = Vn*H1n 
Ja(n,4) = Vn*Vn*H0n 

c wri te(7,*) 'n=' ,n, ' ja(n,4)=', ja(n,4) 

11 Continue 

WRITE(6,*) ' Completed. ' 

WRITE(6,*) ' • 

URITE(6,200) nl ,Hmax 
C 

C Loading the Pink Matrix (size Nl x Mmax) 

C 

DO 22 m=0 # N1 

URITE(6,300) m,cr 
DO 22 k=0,Mmax 

Pmk<m # k)=(0. 0,0.0) 

DO 33 n=0,N2 

S1=( Ja(n, 1 )*T(m,n)*T(k,n))/( I (n,2)*Ja(n,2)) 

Pmk(m,k)=Pmk(m,k)+S1 
33 CONTINUE 

22 CONTINUE 

WRITE<6,*) # Completed. ' 

DO 44 m=0,Nl 

DO 44 k=0,N1 

If (m.EQ.k) Then 

Amk(nt+1 ,k+1 )=Ja(k,3)*Pmk(m,m)- Ja(k,4)*I <m, 1 ) 

ELSE 

Amk(n*1 ,k+1 )=Ja(k,3)*Pmk(m,k) 

ENDIF 



44 


CONTINUE 






100 


FORMAT ( ' 


Loading Ja matrix of size', 13,' x 4 ') 




200 


FORMAT (' 


Loading Pmk Matrix of size', 13,' x',13,' 


.') 


300 


FORMAT ( 1 


Calculating row',i3,a,\) 






RETURN 






C 


END 






c ******************************************************************** 


v- 

r 


SUBROUTINE 


BLOAD(FMHZ,M1,Pmk,Ja,B) 




L 

c 


*********************************************************** 


c 


* 




* 


c 


* 


Variable Definitions 


* 


c 


* 




* 


c 


*********************************************************** 


c 


B(m) = 


Driving Matrix 




c 


Nl 


Number of Coefficient Expansions in Region 1 




c 


I3(m) = 


Trig. Moment Integrations 




c 


FMHZ 


Driving voltage frequency in megahertz 




c 


eps = 


8.854E* 12 




c 


j 


Imaginary One 





C 

Real d, 1(0:500, 1 :3),Pi ,eps,w 
Integer k,M,M1 ,N1 ,Nbig 

Complex B(Nbig), j,Ja(0: 500,1 :4),Pmk(0: 60, 0:500), Sun 
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C0MM0N/Blk1/I , j,Nl ,Pi ,d,Nbig 
eps=8.854E-12 
w=2.0E6*Pi*FMHZ 
DO 22 m=0,Nl 

Sun=(0. 0,0.0) 

DO 33 k=0,Ml 

SLiTFSunH-((I(k,3)\Ja(k,3)*Pmk(m,k))/<I(k,1)*Ja(k,4))) 
33 CONTINUE 

B ( m 1 ) = j *w* eps*Sun/d 
22 CONTINUE 
RETURN 
END 

******************************************************************** 



FUNCTION CSR(R) 

REAL R 
COMPLEX CSR 

Evaluating SQRT(R) with special consideration to neg. reals 
IF(R.GE.O.O) then 
CSR=SQRT(R) 

else 

CSR=(0.,-1.)*SQRT(-R) 

end if 
RETURN 
END 



******************************************************************** 



SUBROUTINE BES1 (Z, JO, J1 ) 

Computing Bessel Functions for n=0,1 with 
Complex Argunent, Z. Direct Power Series Method for 
CABS(Z) .LE. 6 and Hankel's Asymptotic Formula for 
CABS(Z) .GT. 6. Written 11/5/87 by M.A. Morgan 

INTEGER M,M2 

REAL C(34),DM,F(34),G0,P(34),Pi,P2 
COMPLEX Z,Z2,Z3,Z4, JO, J1 , AM,CL,P0,P1 ,Q0,Q1 ,C0,C1 ,S0,S1 
Pi=3. 1415927 
P2=2.0/PI 

I F(CABS(Z) . LE .6.0) THEN 

Utilizing the Direct Power Series Method 

G0= 1.781072 
Z2=0.5*Z 
CL=CLOG(GO*Z2) 

Computing F(m) = m ! and P(m) = 1 ♦ 1/2 ♦ 1/3 ♦ ♦ 1/m 

F ( 1 )=1 .0 
PC 1 )=1 -0 
DO 11 M=2,34 

F(M)=M*F(M-1 ) 

P(M)=*P(M-1)-M.0/M 
11 CONTINUE 

C 

C Computing Power Series Coefficients 

C 

DM=- 1 .0 
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DO 22 M=1,34 

C(M)=DM/(F(M)*F(M)> 

DH=-DM 

22 CONTINUE 

C 

Computing JO and J1 

J0=(1.,0.) 

J1=(0.,0.) 

M=0 
M=H+1 
M2=2*H 

AM=C(M)*(Z2**M2) 

J0=J0+AM 
J1=J1-M*AM 

IF((CABS(AM).GT.1.0E-10>.AND.(M.IT.34)) GO TO 33 
J1=J1/Z2 
return 

ELSE 

Hankel' Asymptotic Formula (Abram. & Stegun p. 364) 

Z2=Z*Z 
Z3=Z*Z2 
Z4=Z*Z3 

P0=1 .0- .0703125/Z2+.1121521/Z4 
Q0=- .125/Z+.0732422/Z3 
P1=1 .0+.1 171 875/Z2- .1441 956/Z4 
Q1=.375/Z- . 10253906/Z3 
C0=CC0S ( Z - . 25*P I ) 

S0=CSIN(Z-.25*PI) 

C1=CCOS(Z-.75*PI) 

S1=CSIN(Z-.75*PI) 

AM=CSQRT(P2/Z) 

J0=AH*(P0*C0-Q0*S0) 

J1=AM*(P1*C1 -Q1*S1 ) 

END IF 
RETURN 
END 



******************************************************************** 



SUBROUTINE HAN1 (Z # H0, Hi ) 

C 

C Computing Hankel Functions for n=0,1 with 

C Complex Argument, Z. Direct Power Series Method for 

C CABS(Z) .LE. 5 and Hankel 's Asymptotic Formula for 

C CABS(Z) .GT. 5. Written 11/6/87 by M.A. Morgan 

C 

INTEGER M,M2 

REAL C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX Z,Z2,Z3,Z4,J0,J1,Y0,Y1,AM,CL,P0,P1,Q0,Q1 

COMPLEX E0,E1,X0,X1,H0,H1,j 

PI=3. 1415927 

P2=2.0/PI 

j=(0. , 1 . ) 

I F ( CABS ( Z ) . LE . 5 . 0 ) THEN 
C 

C Direct Power Series Method 

C 



G0= 1.78072 
Z2=0.5*Z 
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CL=CLOG(G0*Z2) 



Computing F(m) = m I and P(m) = 1 ♦ 1/2 + 1/3 + ....+ 1/m 

F(1 )=1 .0 
P(1)=1.0 
DO 11 M=2,34 

F(M)=M*F(M-1 ) 

P(M)=P(M-1)+1.0/M 

CONTINUE 

Computing Power Series Coefficients 

DM=*1 .0 
DO 22 M=1 ,34 

C(M)=DM/(F(M)*F(H)) 

DM=-DM 

CONTINUE 

Computing JO and J1 



J0=(1.,0.) 

J1=(0.,0.) 

M=0 

M=M+1 

M2=2*M 

AH=C(M)*(Z2**M2) 

J0=J0+AM 

J1=J1-M*AM 

I F((CABS(AM).GT .1 .0E-10).AND.(M.LT .34)) GO TO 33 
J 1 = J 1/Z2 

Computing Y0 and Y1 



H=0 

Y0=CL*J0 

Y1=Z2*CL*J1-0.5*J0 
44 M=H+1 

M2=2*M 

AH=C(M)*P(M)*(Z2**M2) 

Y0=Y0-AM 

Y1=Y1+M*AM 

IF((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 44 

Y0=P2*Y0 

Y1=P2*Y1/Z2 

H0=J0- j*Y0 

H1=J1- j*Y1 

RETURN 

ELSE 

C 

C Hankel' Asymptotic Formula (Abram. & Stegun p. 364 
C 

Z2=Z*Z 

Z3=Z*Z2 

Z4=Z*Z3 

P0=1 .0- .0703125/Z2+. 1 121521/Z4 
00=- . 125/Z+.0732422/Z3 
P1=1 .0+ . 1 171875/Z2- . 1441956/Z4 
01=.375/Z- . 10253906/Z3 
X0=(Z-.25*PI ) 

X1=(Z-.75*PI ) 

E0=CEXP(- j*X0) 

E1=CEXP(- j*X1) 

AH=CSQRT(P2/Z) 
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HO=AH*(PO- j*QO)*EO 
H1=AH*(P1 - j*Q1 )*E1 

ENDIF 

RETURN 

END 



C 

C ******************************************************************** 

c 

SUBROUTINE FACTOR (A,P,N,NMX) 

C PERFORMING ^-DECOMPOSITION WITH PIVOTING ON THE A-MATRIX 
C WRITTEN BY M.A. MORGAN 
C 

INTEGER I / J,JP1,K,P(NMX) > PJ,PR,R,RM1,RP1 

REAL DMAX.ELMAG 

COMPLEX A(NMX,NMX),D<120) 

DO 60 R=1 ,N 

DO 10 K=1 ,N 

10 D(K)=A(K,R) 

RM1=R-1 

IF(RMI.LT.I) GO TO 31 
DO 30 J=1,RM1 
PJ=P(J) 

A(J,R)=D(PJ) 

D(PJ )=D( J ) 

JP1=J+1 
DO 20 I=JP1,N 

20 D(I)=D(I)-A(I,J)*A(J,R) 

30 CONTINUE 

31 CONTINUE 
DMAX=CABS(D(R)) 

P(R)=R 

RP1=R+1 

IF(RPI.GT.N) GO TO 41 
DO 40 I=RP1,N 

ELMAG=CABS(D(I ) ) 

IF(ELMAG.LT.DMAX) GO TO 40 
DMAX=ELMAG 

40 P(R)=I 

41 CONTINUE 
PR=P(R) 

A(R,R)=D(PR) 

D(PR)=D(R) 

IF(RPI.GT.N) GO TO 51 
DO 50 I =RP 1 ,N 

50 A(I,R)=D(I)/A(R,R) 

51 CONTINUE 

60 CONTINUE 

RETURN 

END 

C 

Q ******************************************************************** 

c 

SUBROUTINE SOLVE (A.B.P.N.NMX) 

C BACK SUBSTITUTION TO INVERT THE LINEAR SYSTEM 

C WRITTEN BY M.A. MORGAN 

C 



COMPLEX A ( NMX ,NMX),B(NMX),Y(120), SUM 
INTEGER I , IP1, J,K,P(NMX),PPI 
DO 20 1=1 ,N 
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PP!=P( I ) 

Y( I )=B(PPI ) 

B(PPI )=B( I ) 

IP1=I+1 

IF(IPI.GT.N) GO TO 11 
DO 10 J=IP1,N 

B( J)=B( J)-A(J, I )*YU ) 

10 CONTINUE 

11 CONTINUE 

20 CONTINUE 

DO AO K=1 ,N 

I=N-K+1 

SUM=(0.,0.) 

IP1=I+1 

IF( IP1 .GT.N) GO TO 31 
DO 30 J=IP1,N 

SUM=SUH+A(I,J)*B(J) 



30 


CONTINUE 


31 


CONTINUE 




B(I )=(Y( I )*SUH)/A(I , I ) 


40 


CONTINUE 




RETURN 




END 



******************************************************************* 
PROGRAM PLOT 

Program to Plot a Solid Line Plot with Option to Overlay a 
♦♦++ Line Plot for Comparison. Uses PLOT2 Data File Format. 

Original MS-FORTRAN Version of PLOT2 4/24/87 by M.A. Morgan. 
Hods: 5/30/87 XHIN; 11/3/87 RM-FORTRAN 2nd Plot Option. 

INPUT DATA FORMAT 

Unit #3: Solid Line Plot 

TITLE1- 64 Space Header 
N - # Data Points 15 (Integer*2) 

XMIN - Real Min X value E12.3 

XMAX - Real Max X value E12.3 

F1(N) - Input Data Array E 1 2.3 

C Unit #4: + Symbol Plot 

C 

C TITLE2- 64 Space Header 

C N - # Data Points 15 (Integer*2) 

C XMIN - Real Min X value E12.3 

C XMAX - Real Max X value E12.3 

C F2(N) - Input Data Array E12.3 

C 

CHARACTERS YN, YN1 ,DUM,YN2, SYMBOL 

CHARACTERS LINE 

CHARACTERS CROSS 

CHARACTER‘16 LTIT.CTIT 

CHARACTER*32 SCALE 

CHARACTER*64 TITLE, FNAME 

REAL X( 1025), F1(1 025), F2( 1025) 

INTEGERS N , JROW, JCOL , I SYMBL , I TYPE 
SCALE='is the Vertical Scale Multiplier' 

Line*'--- ' 

CROSS*' +++ ' 
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CALL QBEEP 

URITE(*,*) 'Press RET to Exit from Screen Plot' 

WRITE(*,*) 'Printer Hardcopy ? (Y/N): ' 

READ(*,100) YN 

WRITE(*,*) 'Enter Data File Name for Plot: ' 

READ(*, 100) FRAME 
C Reading from Unit # 3 
0PEN(3,FILE=FNAME) 

READ13, 100) TITLE 
READ(3,120) N 
READ(3, 130) XMIN 
READ13, 130) XMAX 

URITE(*,*) ' Plot TITLE, N, XMIN, XMAX: ' 

WRITE(*,*) TITLE 
WRITEl*,*) N, XMIN, XMAX 

WRITE(*,*) 'Enter Plot Caption (16 char max): ' 

READ(*, 100) LTIT 
DX=(XMAX-XMIN)/(N-1.0) 

FMIN=0.0 

FMAX=0.0 

DO 22 K=1,N 

X(K)=XMIN+(K-1.0)*DX 

READ(3,130) F1(K) 

lF(FKK).LT.FMIN) FMIN=F1(tO 

IF(F1(K).GT.FMAX) FMAX=F1(K) 

22 F2(K)=0.0 

WRITE(*,*) 'Comparison Plot (++++) on Same Graph ??? (Y/N) 
READ(*, 100) YN1 

IF((YN1.EQ.'N').OR.(YN1.EQ.'n')) GO TO 11 
URITE(*,*) 'Enter Data File Name for ++++ Plot: ' 

READ(*, 100) FRAME 
C Reading from Unit # 4 

OPEN ( 4 , F I LE= FRAME ) 

READ(4, 100) TITLE 
READ(4,120) N 
READ(4, 130) XMIN 
READ(4, 130) XMAX 

WRITE(*,*) '++++ Plot TITLE, N, XMIN, XMAX:' 

WRITE(*,*) TITLE 
WRITE(*,*) N, XMIN, XMAX 

WRITE(*,*) 'Enter ++++ Plot Caption (16 char max): ' 

READ(*, 100) CTIT 
DO 27 K=1 ,N 
READ(4, 130) F2(K) 

IF(F2(K).LT.FMIN) FMIN=F2(K) 

IF(F2(K).GT.FMAX) FMAX=F2(K) 

27 CONTINUE 

CLOSE (4) 

11 CONTINUE 

CALL QBEEP 

WRITE(*,*) 'Enter TITLE for Plot (64 Char. Max):' 
WRITE(*,*) ' ' 

READ (*,100) TITLE 
CLOSE (3) 

IF(FMIN.GT.O.O) FMIN=0.0 
IF(FMAX.LT.O.O) FMAX=0.0 
C Computing Scale Factors for Vertical Axis 

ABSMIN=ABS(FMIN) 

A8SMAX=ABS ( FMAX ) 

YMAX=AMAX1 ( ABSM I N , ABSMAX ) 

NSCL=INT(LOG10(YMAX) ) 

IF (YMAX.LT.1.0) NSCL=NSCL-1 

YSCL=10.**NSCL 

FMIN=FMIN/YSCL 
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FHAX=FKAX/YSCL 
ABSM I N=ABSM I N/YSCL 
ABSMAX=ABSMAX/YSCL 
DO 33 K=1 ,N 
F1(K)=F1(K)/YSCL 
33 F2(K)=F2(K)/YSCL 
YMIN=0.0 

IF(FMIN.EQ.O.O) GO TO 37 
35 YMIN=YMIN+1 .0 

IF(ABSMIN.GT.YNIN) GO TO 35 
YM I N=YM I N* FM I N/ABSH I N 
37 CONTINUE 

YHAX=0.0 

IF(FMAX.EQ.O.O) GO TO 41 
39 YHAX=YMAX+1 .0 

IF(ABSMAX.GT.YMAX) GO TO 39 
YHAX=YMAX* FKAX/ABSHAX 
41 CONTINUE 

C Calling GRAFKATIC Routines and Plotting FI Solid Line Graph 
ITYPE=1 
ISYHBL=-2 
CALL OSHODE(6) 

CALL OCMOV(0,21) 

WRITE(*,150) YSCL 

CALL QPTXT(32, SCALE, 3,25, 20) 

CALL QPLOT(100,530, 12, 147,XMIN,XMAX,YMIN,YHAX,0. ,0. ,0, 1 .,1.5) 
CALL QSETUP(0,3, ISYMBL.3) 

XMA JOR=XMAX/5 . 0 

CALL QXAXIS(XMIN,XMAX,XMAJOR, 1 ,1,2) 

YMAJOR=1 .0 

CALL QYAXIS<YMIN,YMAX,YMAJOR, 1 ,1,2) 

CALL QPTXT(64, TITLE, 3, 14,24) 

CALL QPTXT(4,LINE,3,14,22) 

CALL QPTXT(16,LTIT,3,18,22) 

JROW=(ABS(YMIN)/(ABS(YMAX)+ABS(YMIN)))*135 

JCOL=70-430*XMIN/(XMAX-XMIN) 

CALL QGTXT13, '0.0' ,3, JCOL, JR0W,0) 

CALL QTABL(ITYPE,N,X,F1) 

IF((YN1.EQ.'N').OR.(YN1.EO.'n')) GO TO 43 
C Plotting F2 Graph using + Symbol 

CALL QPTXT(7,CROSS,3,34,22) 

CALL QPTXT(16,CTIT,3,41 ,22) 

ISYMBL=43 

ITYPE=0 

CALL OSETUP(0,3, ISYHBL,3) 

CALL QTABL(ITYPE,N,X,F2) 

43 CONTINUE 
READ(*, 100) DUH 

I F( (YN.EQ. 'N' ) .OR. (YN.EO. 'n' )) GO TO 44 
CALL QPSCRN 

44 CONTINUE 

CALL QSMOOE(2) 

WRITEC*,*) 'End of PLOT Program' 

100 FORMAT(A) 

120 FORHAT(I5) 

130 F0RMAT(E12.3) 

150 FORMAT(13X, 1PE1 1 .3) 

STOP 

END 
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APPENDIX C 



SOFTWARE FLOW CHART FOR MONO WITHOUT GAP FIELD EXPANSION 
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ALOAD 

Subroutine 
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APPENDIX D 



SOURCE CODE FOR MONO WITHOUT GAP FIELD EXPANSION 

PROGRAM MONO 
C (Revision 8.1) 

C 

C Mod of Mono 7.2, not using the expansion for the Gap Region 
C Computing Currents on a Monopole Between Parallel 
C Plates using Cylindrical Harmonic Expansions in 2 Regions. 

C Links with SUBS4.08J . Created 05/88 LT. R.C. Hurley 



C Modified 10/88 M. A. Morgan 

C 

C Variables 

C FMHZ = Frequency in MHz 

C h = Monopole height in Meters 

C a = Monopole radius in Meters 

C d = Distance between lower ground plane and feed point 

C in Meters 

C l = Distance between ground planes in Meters 

C q = height from top of monopole to upper ground plane 

C in Meters 

C j = Imaginary One 

C Pi = Value of pi (3.1415927) 

C K0 = 2*pi*F/c 

C where F = frequency of operation in Hz 

C c = Speed of Light 

C K2 = K0**2 

C N1,N2 = Nunber of Coefficients in each of the two regions 

C NBIG =• Size of the arrays used in FACTOR and SOLVE 

C Unit 0 = Standard Input Device 

C Unit 6 = Screen Output 

C Unit 7 = Main Output File to Disk (Name selected by User.) 

C ***** Variable Definitions ***** 

C 

$Large:Amk 



INTEGER k,m,mpn,n,Nl ,N2,Nbig,NI , Nmax,Ns,NZ,P(61 ) 

REAL a,ai , Am, An, c, Cl ,CM,CP,CR,d,di ,Del taZ,Denom,Dz,EzM,eps, FMHZ 

REAL h,hi , 1(0:60, 1 :3),K0,K2, l , l i ,Lambda,Pi ,q,Rin,S 

REAL T (0:60,0: 60 ),w,Xin,y,yhat,Z,Zmin,Zmax 

COMPLEX Amk(61 ,61 ) ,B(61 ) ,Cl ,C2, IZ(61 ) , j , JO, J1 , Ja(0:60, 1 :4) 

COMPLEX Pmk(0 :60, 0:60), Sum, Y0,Y1,Zin 

CHARACTERS Fname 

CHARACTERS Bell 

CHARACTER*! HOW 

CHARACTER*! INPUT 

CHARACTER*! Eok 

CHARACTER*12 INname 

COMMON/Blkl/I , j ,N1 ,Pi ,d,Nbig 

COMMON /B l k2/K2, l , Nmax,N2,q,a,T 

C 

c ********************************* 

C * Main Program Code * 

r ********************************* 
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noon 



C 



* Input Section * 
*********************** 



Bel l=CHAR(7) 

Nbig=61 

c=3.0E+08 

WRITE(6,*) 'Do you wish to use an input data file [Y/N]?' 

READ(0, 100) HOW 

IF (HOU.EQ.'y'.OR.HOW.EQ.'Y') THEN 

WRITE(6,*) 'Enter input data file name with extension 
READ(0, 100) INname 
OPEN (Uni t=4 , F i l e= I Nname ) 

READ(4, 100) INPUT 

IF ((INPUT .EQ. 'w') .OR. (INPUT .EQ. 'U') .OR. (INPUT .EQ. 

+ 'm') .OR. (INPUT .EQ. 'M')) THEN 

READ(4,*) FMHZ,hi,ai,di,li,Nl,N2,NI 

ELSE 

WRITE(*,*) Bell 

WRITE(6,*)'Error detected in first line of data file ' 

+, INname 

WRITE(6,*) 'First line must be either M or W for meter 
+ values or wavelength values.' 

STOP 

END IF 

ELSE 

2 URITE(6,*) 'Enter data in either [M]eters or [U]avelength [M 

+/U]?' 

READ(OJOO) INPUT 

IF (INPUT. EQ.'m'. OR. INPUT. EQ.'M') THEN 

WRITE(6,*) 'Enter Frequency in MHz: ' 

READ(0,*) FMHZ 

URITE(6,*) 'Enter Monopole Height, h, in meters: ' 

READ(0,*) hi 

WRITE(6,*) 'Enter Monopole Radius, in meters: ' 

READ(0,*) ai 

WRITE(6,*) 'Enter Feed Gap, d, in meters: ' 

READ(0,*) di 

WRITE(6,*) 'Enter Upper Plane Height, l (meters): ' 

READ(0,*) li 

ELSEIF (INPUT. EQ.'w'. OR. INPUT. EQ.'W') THEN 

WRITE(6,*) 'Enter Frequency in MHz: ' 

READ(0,*) FMHZ 

WRITE(6,*) 'Enter Monopole Height, h, as factor of wav 

+elength : ' 

READ(0,*) hi 

WRITE(6,*) 'Enter Monopole Radius, as factor of wavele 

+ngth : ' 

READ(0,*) ai 

WRITE(6,*) 'Enter Feed Gap, d, as factor of wavelength 

+ : ' 

READ(0,*) di 

WRITE(6,*) 'Enter Up>per Plane Height, l, as factor of 
+wavelength : ' 

READ(0,*) li 

ELSE 

WRITE(*,*) Bell 

WRITE(6,*) 'You must enter either an M or W. 

+Reenter.. ' 

GOTO 2 

END IF 

URITE(6,*) 'Enter # Region 1 Coeffs. (N1 .LE. 60): ' 

READ(0,*) N1 
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WRITE(6,*) 'Enter # Region 2 Coeffs. (N2 .IE. 60): ' 

READ(0,*) N2 

WRITE(6,*) 'Enter No. Points for 1(2) including ends: ' 

READ(0,*) NI 

ENDIF 



C *********************************************** 

* * 

* Calculate the Physical Parameters * 

* in both Wavelenght and Meters * 

* * 

C *********************************************** 



Lambda=c/(FMHZ*10**(6)) 

IF (INPUT. EQ.'W'. OR. INPUT. EG.'W') THEN 
h=hi*Lambda 
a=ai*Lambda 
d=di*Lambda 
1=1 i*L ambda 

ELSE 

h=hi 

hi=h/ Lambda 
a=ai 

ai = a/ Lambda 
d=di 

di =d/ Lambda 
1 = 1 i 

l i = l/Lambda 

ENDIF 

C *********************************************** 

C * * 

C * Is Upper Plate Height an * 

C * Even Multiple of Wavelength? * 

C * * 

Q *********************************************** 



IF (MOO (( 1*2), Lambda ).LT. IE-05) THEN 
WRITE(*,*) Bell 
WRITE(6,*) ' ' 

WRITE(6,*) 'ERROR! Upper Plate Height is almost an even mul 
+tiple of half the wavelength.' 

WRITE(6,*) 'Ensure l is NOT within 0.00001 of an integer mul 
♦tiple of half the wavelength.' 

STOP 

ENDIF 



C 

C 

C 

C 

C 

C 

C 



*********************************************** 
* * 

* Is Upper Plate Height an * 

* Even Multiple of Antenna Height? * 

* * 
*********************************************** 



IF ((MOD(l,h)).LT. IE-05) THEN 
WRITE(*,*) Bell 
WRITE(6,*) ' ' 

WR1TE(6,*) 'ERROR! Upper Plate Height is almost an even mul 
+tiple of the antenna height.' 

WRITE(6,*) 'Ensure l is NOT within 0.00001 of an integer mul 
♦tiple of h (in wavelength).' 

STOP 

ENDIF 

WRITE(6,*) 'Enter output file name and extension: ' 
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READ(O.IOO) Fname 
OPEM (Uni t=7, Fi le=Fname) 

WRITE(7,210) FMHZ 
WRITE(7,220) h,hi 
WRITE(7,230) a,ai 
URITE(7,240) d,di 
WRITE(7,250) l, If 

WRITE(7,*) 'Region I Coeffs. (Ml): \N1 
WRITE(7,*) 'Region II Coeffs. (N2): ',N2 
WRITE(7,*) 'No. Points for 1(2): ',NI 

j=(0. 0,1.0) 

Pi=3. 1415927 

K0=Pi*FMHZ/150.0 

K2=K0*K0 

eps=8.854E-12 

w=2.0E6*Pi*FMHZ 

yhat=w*eps 

q=l-h 

Nmax=MAX0(Nl ,N2) 

WRITE(6,*) 'Loading the Moment Matrices ' 

WRITE(6, *) ' ' 

Q *********************************************** 

c * * 

C * Evaluating Moment Integrations * 

C * Over Region h < z < l. * 

C * * 

Q *********************************************** 
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C 

C 

C 

C 

C 

C 



T(0,0)=q 
DO 11 k=1,N2 

T ( 0 , k )=0 . 0 
CONTINUE 
DO 22 m=1,Nmax 
Am=m*Pi/l 

DO 22 r>=0,N2 

Ar>=n*Pi/q 

IF((ABS(Am-An)).GT.1.0E-3) THEN 
D enom=Am*Am- An* An 
S=SIN(Am*h) 
T(m,n)=-Am*S/Denom 

ELSE 

mpn=m+n 

T(m, n)=0.5*q*((-1)**mpn) 

END IF 
CONTINUE 

WRITE(6,*) ' Moment Matrix T Loaded.' 

URITE(6,*) ' ' 



*********************************************** 



* Loading the Im Matrix for the * 

* Largest Number of Coefficients * 

* * 

*********************************************** 



DO 33 m=0,Nmax 

IF(m.EQ.O) Then 
I(m,1)=l 
I (m,2)=q 
I(m,3)=d 

ELSE 

I(m,1)=l/2.0 

I(m,2)=q/2.0 
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Am=m*P i / l 

I(m,3)=SIN(Am*d)/Am 

END IF 

33 CONTINUE 

WRITEC6,*) ' Moment Matrix I Loaded.' 

WRITE<6,*) ' ' 

WRITE(6,*) 'Loading System Matrix ' 

CALL ALOAD(Ja,Pmk,Amk) 

WRITE(6,*) 'System Matrix Amk Loaded.' 

WRITE(6,*) ' ' 

WRITE(6,*) 'Loading Driving Vector ' 

DO 40 m=0,N1 

B(m+1 )=j*yhat*I (m,3)/d 
40 CONTINUE 

WRITE(6,*) 'Driving Vector B Loaded.' 

WRITE(6,*> ' ' 

WRITE(6,*) 'Solving the Linear System .. Please Wait .... ' 
Ns=N1+1 

CALL FACTOR ( Amk, P, Ns, Nbig) 

CALL SOLVE(Amk,B,P,Ns,Nbig) 

WRITE(7,*) 'Listing the Expansion Coefficients' 

WRITE(6,*) ' Listing the Expansion Coefficients ' 

WRITE(7,*) 'N and B(N): ' 

DO 44 N=0,N1 

WRITE (7,*) N,B(N+1) 

44 CONTINUE 

*********************************************** 

★ * 

* Computing Current from the * 

C * Expansion Coefficients * 

C * * 

C *********************************************** 

WRITEC6,*) ' # 

WRITE(6,*) 'Computing current on the monopole ' 

WRITE(7, *) 'Listing the Current on the Monopole' 
DeltaZ=(h-d)/(NI -1 ) 

WRITE(7,*) 'Position Distance 1(2) (A) 

WRITEC7,*) ' Nunber (Meters) Mag 

y=d 

k=1 

55 $un=(0.,0.) 

DO 66 n=0,N1 

Sun=Sum+(B(n+1)*Ja(n,3)*cos(n*Pi*y/l)) 

66 CONTINUE 

IZ(k)=2.0*Pi*a*$un 
CM=CABS( IZ(k)) 

CR=REAL(IZ(k)) 

CI=AIMAG( IZ(k)) 

CP=180.*ATAN2(CI,CR)/pi 
WRITE(7 f 200) k,y,CM,CP 
y=y+DeltaZ 
k=k+1 

IF(y.LE.(h+1.E-6)) GO TO 55 
WRITE(6,*) 'Complete.' 

Q *********************************************** 

C * * 

C * Computing Input Impedance * 

c * * 

c *********************************************** 



Zin=(1. 0,0.0)712(1) 



Phase' 
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c 

c 

c 

c 

c 

c 

c 

c 

c 



88 



99 



C 

C 

C 

C 

C 

C 



Rin=REAL(Zf n) 

Xin=AIMAG(Zin) 

WRITE(7,*) 'Input Resistance = ',Rin 
WRITE(7,*) 'Input Reactance = ',Xin 

*********************************************** 
★ * 

* Calculate Electric Field * 

* Distribution * 

* * 
*********************************************** 



Checking Ez(a,z) of solution at 200 points from 0 to L 
WRITE(6,*) 'Do you want to compute the Electric Field (Y/N)? ' 

READ(0, 100) Eok 

IF ((Eok .EQ. 'Y') .OR. (Eok .EQ. 'y')) THEN 
NZ=200 
ZMIN=0.0 
ZMAX=L 

DZ=ZMAX/(NZ- 1 .0) 

OPEN(3,FILE='EZ.DAT') 

WRITE(3, 100) 'EZ vs. Z' 

URITE(3, 102) NZ 
WRITE(3, 103) ZHIN 
URITE(3, 103) ZMAX 
DO 99 K=1,NZ 

Z=(K-1.)*DZ 

C1=(0.,0.) 

DO 88 N=0,N1 

C1=C1+B(N+1)*Ja(N,4)*COS(N*Pi*z/l) 

CONTINUE 

C1=C1/(j*yhat) 

EzM=CABS(Cl ) 

WRITE(3, 103) EzM 

CONTINUE 

WRITE(6,*) ' Electric Field data stored in file EZ.DAT. ' 

ELSE 

WRITE(6,*) ' Completed.' 

END IF 

*********************************************** 



* I/O Format Statements * 

* * 
*********************************************** 



100 


FORMAT (A) 








102 


F0RHATO5) 








103 


FORMAT(F12.3) 








200 


FORMAT(2X, I3,5X,F8.3,5X,1PE12.3,5X,0PF8.2,1PE12.3,5X f 0PF8.2, 




+1PE12.3,5X # 0PF8.2) 








210 


FORMAT (' Driving Frequency, 


' ,9X, F7.3, ' MHz') 






220 


FORMAT (' Monopole Height, 


' ,9X, F7.3, ' meters 


( ' f F5.3, ' 


Lambda)' ) 


230 


FORMAT ( ' Monopole Radius, 


' ,9X, F7.3, ' meters 


('.F5.3,' 


Lambda)' ) 


240 


FORMAT ( ' Feed Gap Distance, 


',9X,F7.3,' meters 


('.F5.3,' 


Lambda)' ) 


250 


FORMAT ( ' Upper Plane Height, 
STOP 


',9X,F7.3,' meters 


('«F5.3,' 


Lambda)' ) 




END 
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Subroutine Sorce Code for MONO without Gap Field Expansion 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



SUBROUTINE ALOAD( Ja,Pmk,Amk) 

(Revision 8.1) 

******************************************** 

* * 

* Variable Definitions * 

* * 

******************************************** 

l = Distance between ground planes 

Radius = Radius of monopole 

q = Distance between top of monopole and upper plane 
d = Antenna feed gap 

K2 = K0**2 (From main program) 

T(m,n) = Moment integration result as function of m and n 

I(m # x) = Moment integration result of order x given m 

j = Imaginary One 

z = Intermediate value for Bessel Subroutine 

un, vn = Eigenvalues 

JO = Bessel Function of the First Kind of Order 0 

J1 = Bessel Function of the First Kind of Order 1 

HO = Hankel Function of the Second Kind of Order 0 

Hi = Hankel Function of the Second Kind of Order 1 

SUM = Summation over n for the Amk matrix 

SI = Individual term used in SUM 

m = Colum position of matrix Amk 

k = Row position of matrix Amk 

N1 = Nunber of subdivisons in Region 1 

N2 = Nimber of subdivisons in Region 2 

Amk = Complex Matrix value that describes the system 

This SUBROUTINE loads the complex A matrix that describes the system 
of size N1 by Nl. Maximum size of matrix is set at 50 by 50. This 
SUBROUTINE calles to the following subroutines found in SUBS. FOR: 

CSR - Complex square root that preserves negative imaginary 
portion. 

BES1 - Computes the Bessel Function values for a given z. 
Tmint - Computes the value of the moment integratation T 
Imint - Computes the value of the moment integratation I 
Written 5/88 by R.C. Hurley using Microsoft Fortran 4.01. 



$LARGE:Amk 



INTEGER k,m,n,Nbig,Nmax,N1 ,N2 
REAL l ,a,q,d,K2,T(0:60,0:60), I (0:60, 1 :3),Pi ,R1 ,R2 
COMPLEX j,zv, zu,vn,un, JOn, Jin, H0n, H1n,SUM,Sl 
COMPLEX CSR, Amk(Nbig,Nbig),Ja( 0:60, 1:4) 

COMPLEX Pmk(0:60,0:60) 

CHARACTER*! cr 

COMMON/B lkl/1 , j ,N1 ,Pi ,d,Nbig 
COMMON/B l k2/K2 , l , Nmax , N2 , q, a , T 
cr=CHAR( 13) 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



This section loads a matrix consisting of Bessel fuctions and 
other values required to load the Amk matrix. The colums are 
f i l led as follows: 

Ja(n,1) = un*Jo(un*a) in numerator of Sum 

Ja(n,2) = Jl(un*a) in denominator of Sum 

Ja(n,3) = vn*Hl(vn*a) used for all locations 

Ja(n,4) = vn*vn*Ho(vn*a) used in diagonal elements 
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c 

c 

c 



33 

22 



44 

100 

200 

300 



WRITE(6, 100) Nmax 
DO 11 n=0,Nmax 

R1=K2-((n*Pi/q)*(n*Pi/q)) 

R2=K2-((n*Pi/l)*(n*Pi/l )) 

un=CSR(R1) 

vn=CSR(R2) 

zu=un*a 

zv=vn*a 

Call BESKzu, J0n,J1n) 

Call HAN1(zv,H0n,H1n) 

Ja(n, 1) = Un*J0n 
Ja(n,2) = Jin 
Ja(n,3) = Vn*H1n 
Ja(n,4) = Vn*Vn*H0n 
Cont i nue 

VRITE(6,*) ' Completed.' 

WRITE(6,*) ' ' 

WRITE(6,200) Nl,Nmax 

Loading the Pmk Matrix (size N1 x Nmax) 

DO 22 m=0,N1 

URITE<6,300) m,cr 
DO 22 k=0,Nmax 

Pmk(m,k)=(0. 0,0.0) 

DO 33 n=0,N2 

Sl=(Ja(n, 1 )*T(m,n)*T(k,n))/( I (n,2)*Ja(n,2)) 

Pmk(m, k)=Pmk(m, k)+S1 
CONTINUE 

CONTINUE 

WRITEC6,*) ' Completed. ' 

DO 44 m=0,N1 

DO 44 k=0,Nl 

If(m.EO.k) Then 

Amk(m*1 ,k+1)=Ja(k,3)*Pmk(m,m)-Ja(k,4)*I(m, 1 ) 

ELSE 

AmkCnrH-l ,k+1 )=Ja(k,3)*Pmk(m,k) 

END IF 

CONTINUE 

FORMAT ( ' Loading Ja matrix of size', 13,' x 4 ') 

FORMATS Loading Pmk Matrix of size', 13,' x',13,' ') 

FORMAT ( ' Calculating row',i3,a,\) 

RETURN 

END 



C 

C ******************************************************************** 

C 



FUNCTION CSR(R) 

REAL R 
COMPLEX CSR 

C Evaluating SORT(R) with special consideration to neg. reals 
IF(R.GE.O.O) THEN 
CSR=SQRT(R) 

ELSE 

CSR=(0.,-1.)*SQRT(-R) 

END IF 
RETURN 
END 
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******************************************************************** 



SUBROUTINE BES1(Z, JO, Jl) 

Computing Bessel Functions for n=0,1 with 
Complex Argument, Z. Direct Power Series Method for 
CABS(Z) .LE. 6 and Hankel's Asymptotic Formula for 
CABS(Z) .GT. 6. Written 11/5/87 by M.A. Morgan 

INTEGER M,M2 

REAL C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX Z,Z2,Z3,Z4, JO, Jl ,AM,CL,P0,P1 ,00,01 ,C0,C1 ,S0,S1 
Pi=3. 1415927 
P2=2.0/PI 

IF(CABS(Z).LE.6.0) THEN 

Utilizing the Direct Power Series Method 

G0= 1.781072 
Z2=0.5*Z 
CL=CLOG(GO*Z2) 

Computing F(m) = m ! and P(m) = 1 + 1/2 + 1/3 + ....+ 1/m 

F(1)=1.0 
P(1)=1.0 
DO 11 M=2,34 

F(M)=M*F(M-1) 

P(M)=P(M-1)+1.0/M 
1 CONTINUE 

Computing Power Series Coefficients 



DM=-1 .0 
DO 22 M=1 ,34 

C(M)=DM/(F(M)*F(M)) 

DM=-DM 

CONTINUE 

Computing JO and Jl 



J0=(1.,0.) 

J1=(0.,0.) 

M=0 

33 M=M+1 

M2=2*M 

AH=C(M)*(Z2**M2) 

J0=J0+AM 

J1=J1-M*AM 

IF((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 33 

J1=J1/Z2 

return 

ELSE 

C 

C Hankel' Asymptotic Formula (Abram. & Stegun p. 364) 

C 

Z2=Z*Z 

Z3=Z*Z2 

Z4=Z*Z3 

P0=1 .0- .0703125/Z2+ . 1 121521/Z4 
00=- . 125/Z+.0732422/Z3 
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P1=1 .0+. 1171 875/Z2-. 1441 956/Z4 
Q1=.375/Z- .10253906/Z3 
CO=CCOS(Z-.25*PI> 
$0=CSIN(Z-.25*PI ) 
C1=CC0S(Z-.75*PI) 
$1=C$IN(Z-.75*PI ) 
AM=CSQRT(P2/Z) 
J0=AM*(P0*C0-Q0*$0) 
J1=AM*<P1*C1-01*$1) 

END IF 

RETURN 

END 



******************************************************************** 



SUBROUTINE HAN1 (Z, HO, Hi ) 

Computing Hankel Functions for n=0,1 with 
Complex Argument, Z. Direct Power Series Method for 
CABS(Z) .LE. 5 and Hankel's Asymptotic Formula for 
CABS(Z) .GT. 5. Written 11/6/87 by M.A. Morgan 

INTEGER M,M2 

REAL C(34),DM, F(34),G0,P(34),Pi ,P2 

COMPLEX Z,Z2,Z3,Z4,J0,J1,Y0,Y1,AM,CL,P0,P1,Q0,Q1 

COMPLEX E0,E1 ,X0,X1 ,H0,H1 , j 

PI=3. 1415927 

P2=2.0/PI 

j=(0.,1.) 

I F(CABS(Z) .LE .5 .0) THEN 

Direct Power Series Method 

G0= 1.78072 

Z2=0.5*Z 

CL=CLOG(GO*Z2) 

Computing F(m) = m ! and P(m) = 1 + 1/2 + 1/3 + + 1/m 

F( 1 )=1 .0 
P(1)=1.0 
DO 11 M=2,34 

F(M)=M*F<M-1) 

P(M)=P(M-1)+1.0/M 
1 CONTINUE 

Computing Power Series Coefficients 



DM=-1 .0 
DO 22 M=1 ,34 

C(M)=DM/(F(M)*F(M)) 

DM=-DM 

CONTINUE 

Computing JO and Jl 



J0=(1.,0.) 

J1=(0.,0.) 

M=0 

33 M=M+1 

M2=2*M 

AM=C(M)*(Z2**M2) 
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J0=J0+AM 
J1=J1 -M*AM 

I F((CABS(AM) .GT .1 .0E- 10).AND. (H.LT.34)) GO TO 33 
J1=J1/Z2 
C 

Computing YO and Y1 



M=0 

Y0=CL*J0 

Y1=Z2*CL*J1-0.5*J0 

M=M+1 

M2=2*M 

AM=C(M)*P(M)*(Z2**M2) 

Y0=Y0-AH 

Yl*Y1+M*AM 

IF((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 44 

Y0=P2*Y0 

Y1=P2*Y1/Z2 

H0=J0-j*Y0 

HI = J 1 - j*Y1 

RETURN 

ELSE 

Hankel' Asymptotic Formula (Abram. & Stegun p. 364 

Z2=Z*Z 

Z3=Z*Z2 

Z4=Z*Z3 

P0=1 .0- .07031 25/Z2+. 11 21521/Z4 
Q0=- .125/Z+.0732422/Z3 
P1=1 .0+.1 171875/Z2- .1441956/Z4 
Q1=.375/Z- .10253906/Z3 
X0=(Z-.25*PI) 

X1=(Z-.75*PI) 

E0=CEXP(- j*X0) 

E1=CEXP(- j*X1 ) 

AM=CSQRT (P2/Z) 

H0=AM*(P0- j*Q0)*E0 
H1=AM*(P1 - j*01)*E1 

ENDIF 

RETURN 

END 



******************************************************************** 



SUBROUTINE FACTOR (A,P,N,NMX) 

C PERFORMING LU-DECOHPOSITION WITH PIVOTING ON THE A-HATRIX 
C 

INTEGER I , J, JP1,K,P(NMX),PJ,PR,R,RM1,RP1 

REAL DMAX.ELKAG 

COMPLEX A(NMX,NMX) ,0(120) 

DO 60 R=1,N 

DO 10 K*1,N 

10 D(K)=A(K,R) 

RM1=R-1 

IF(RMI.LT.I) GO TO 31 
DO 30 J*1,RM1 
PJ=P(J) 

A(J,R)=D(PJ) 

D(PJ)=D(J) 

JP1=J+1 
DO 20 I = JPl ,N 
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D(I )=D(I )-A(I , J)*A(J,R) 



20 

30 CONTINUE 

31 CONTINUE 
DMAX=CABS(D(R)) 

P(R)=R 

RP1=R+1 

IF(RPI.GT.N) GO TO 41 
DO 40 I=RP1 ,N 

ELMAG=CABS(D(I)) 

IF(ELMAG.LT.DKAX) GO TO 40 
DMAX=ELMAG 

40 P(R)=I 

41 CONTINUE 
PR=P(R) 

A(R,R)=D(PR) 

D(PR)=D(R) 

IF(RPI.GT.N) GO TO 51 
DO 50 I=RP1,N 

50 A(l,R)=D(l)/A(R,R) 

51 CONTINUE 

60 CONTINUE 

RETURN 

END 

C 

g ******************************************************************** 

c 

SUBROUTINE SOLVE <A,B,P,N,NMX) 

C BACK SUBSTITUTION TO INVERT THE LINEAR SYSTEM 

COMPLEX A(NMX,NMX),B(NMX),Y(120),SUM 
INTEGER I , IP1 , J,K,P(NMX),PPI 
DO 20 1=1 ,N 

PPI=P(I) 

Y( I )=B(PPI ) 

B(PPI)=B(I) 

IP1=I+1 

IF(IPI.GT.N) GO TO 11 
DO 10 J=IP1,N 

B(J)=B(J)-A(J,I)*Y(I) 

10 CONTINUE 

11 CONTINUE 

20 CONTINUE 

DO 40 K=1,N 

I=N-K+1 

SUM=(0.,0.) 

IP1=I+1 

IF(IPI.GT.N) GO TO 31 
DO 30 J=IP1,N 

SUM=SUM+A( I , J)*B( J) 



30 


CONTINUE 


31 


CONTINUE 




B(I)=(Y(I )-SUM)/A( 1,1) 


40 


CONTINUE 




RETURN 




END 
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APPENDIX E 



SAMPLE INPUT/OUTPUT DATA FILES USED WITH MONO 

This Appendix lists the input and output data file 
structures used within MONO. There are two different types 
of input data files, one that is entered at the keyboard at 
run time and one which allows unattended processing of a 
number of runs without user intervention. In each case, the 
data fields may be either arranged in column format with one 
entry per line, or all on one line separated by commas. 
Additionally, a variable type is set by the program which 
allows a real variable to read both an input of "1" or "1." 
as the same. 

A. INPUT DATA FILE STRUCTURE 

This data file may be adjusted to allow the user the 
option of selecting either the output file name or, for the 
case of M0N08, the calculation of the electric field along the 
antenna at run time. Below is a commented data set that could 
be used. Note the comments are NOT to be entered in the data 
set as errors would occur. 
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Data File 
Entry 


Comment 


w 


This identifies the basis for all physical 
measurements of the antenna. Using a "w" 
or "W" signifies all measurements are in 
terms of wavelength, while a "m" or "M" 
indicates measurements are in terms of 
meters. 


299.8 


Driving source frequency in MHz 


.24 


Antenna height, h . Can be entered as 
"0.24". 


.01 


Antenna radius, a. 


.02 


Gap distance, q. 


2 . 17 


Ground plate spacing, /. 


30 


Number of modes for Region I, Nj. 


30 


Number of modes for Region II, N 2 . 


60 


Number of modes for Gap expansion, M x . 
Not used for M0N08. 


30 


Number of points to calculate the current 
along the antenna. 


d:\file.out 


DOS path and name to assign to output 
file. Maximum length is 25 characters. 
This entry may also be entered at the key- 
board at run time. 


"y" or "n" 


Indicates whether to calculate the elec- 
tric field along the antenna. Test is NOT 
case sensitive. If yes, data is stored as 
1 EZ . DAT' in the current directory in for- 
mat for PLOT. FOR found in Appendix B. 

This may be entered via the keyboard at 
run time. 



The same data file can be arranged as follows: 
w, 299 . 8 , .24, . 0 1. 02, 2. 17, 3 0, 3 0, 60, 30, d:\file. out, n 
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The data file is easily modified to allow the user to 
group a series of runs together to be batch processed later 
using the standard DOS redirection command. For example, to 
use input data file named 'RUNl' and to have the normal screen 
output redirected to a running log file named 'OUTPUT.LOG', 
the correct syntax would be: 

MONO <|>a//i]RUNl »[ path ] OUTPUT . LOG 
where path is the DOS path where the individual files are 
stored. This parameter is optional and defaults to the 
current directory. A number of runs could then be completed 
without user input by combining a series of these commands in 
one batch file. The file structure is the same as that listed 
above except that an ' n' is inserted at the beginning of the 
file to indicate that a separate input file is not to be used. 
Therefore, the same data file used above would be modified as 
seen below: 

n, w, 299 . 8 , .24, .01. 02, 2. 17, 30 , 30 , 60 , 30 , d: \f ile . out , n 

B. SAMPLE OUTPUT FILE 

The output file generated by M0N07 for the inputs used 
above is listed below: (* indicates not included in output 

from M0N08) 
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299.800 MHz 
.240 meters 
.010 meters 
.010 meters 
2.171 meters 
30 
30 
60 
30 



Driving Frequency, 

Monopole Height, 

Monopole Radius, 

Feed Gap Distance, 

Upper Plane Height, 

Region I Coeffs. (Nl) : 

Region II Coeffs. (N2) : 

No. Gap Coeffs. (Ml) : 

No. Points for I(z): 

Listing the Expansion Coefficients 
N and B(N) : 

( -9 . 259474E-05 , - 3 . 712934E- 
(-1 . 716743E-04 , -7 . 285348E- 
(-1 . 250944E-04 , -6 . 843712E- 
(-8 . 932992E-06 , -5 . 970440E- 
( 6 . 372479E-04 , - 3 . 034107E- 
(-7.747178E-04,-5.976821E- 
(-4 . 176759E-04 , -5 . 10618 1E- 
(-2 . 877242E-04 , -4 . 189367E- 
(-2 . 029108E-04 , -3 . 268855E- 
(-1 . 371401E-04 , ~2 . 384746E- 
(-8. 359972E-05,-1.573042E- 
(-4 . 017 067E-05 , -8 . 621241E- 
(-6 . 131331E-06 , -2 . 729695E- 
( 1 . 893293E-05 , 1 . 824955E- 
( 3 . 552895E-05 , 5 . 013411E- 
( 4 . 440288E-05 , 6 . 891159E- 
(4 . 656793E-05, 7 . 587663E- 
(4.326494E-05,7. 290936E- 
( 3 . 588156E-05 , 6 . 228465E- 
( 2 . 58803 IE-05 , 4 . 64 9168 E- 
( 1 . 465865E-05 , 2 . 799412E- 
(3 . 497886E-06, 9 . 087322E- 
(-6 . 533222E-06 , -8 . 273393E- 
(-1 . 463069E-05, -2 . 258123E- 
(-2 . 028524E-05 , -3 . 28428 2E- 
(-2 . 328759E-05 , -3 . 859688E- 
(-2.370753E-05,-3. 9882 46E- 
(-2 . 185403E-05 , -3 . 717663E- 
(-1 . 821469E-05 , -3 . 127903E- 
( -1 . 339421E-05 , -2 . 322528E- 
(-8 . 037723E-06 , -1 . 412969E- 
Listing the Current on the Monopole 



( .240 Lambda ) 
( .010 Lambda) 
( .010 Lambda) 
(2.170 Lambda) 



0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 



04) 

04) 

04) 

04) 

04) 

04) 

04) 

04) 

04) 

04) 

04) 

05) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 

05) 

06) 
06) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 



Position 

Number 



Distance 

(Meters) 



Mag 



I(z) (A) 



Phase 



1 


.010 


2.316E-02 


-15.04 


2 


.018 


2 . 309E-02 


-15.98 


3 


.026 


2 . 299E-02 


-17.24 


4 


.034 


2 . 284E-02 


-18.58 


5 


.042 


2 . 262E-02 


-19.80 


6 


.050 


2 . 234E-02 


-20.78 
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7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 



058 


2.200E-02 


-21.49 


066 


2 . 161E-02 


-21.97 


073 


2.121E-02 


-22.34 


081 


2 . 082E-02 


-22.72 


089 


2 . 044E-02 


-23.19 


097 


2 . 007E-02 


-23.78 


105 


1 . 970E-02 


-24.44 


113 


1.929E-02 


-25.11 


121 


1 . 881E-02 


-25.70 


129 


1.822E-02 


-26.16 


137 


1.752E-02 


-26.50 


145 


1 . 669E-02 


-26.74 


153 


1 . 574E-02 


-26.95 


161 


1 . 467E-02 


-27.19 


169 


1.349E-02 


-27.49 


177 


1 . 222E-02 


-27.86 


185 


1 . 088E-02 


-28.27 


193 


9 . 485E-03 


-28.64 


200 


8 . 069E-03 


-28.91 


208 


6.671E-03 


-29.03 


216 


5.33 IE-03 


-28.99 


224 


4 . 091E-03 


-28.84 


232 


2 . 987E-03 


-28.72 


240 


2 . 042E-03 


-28.84 


ice = 


41.703650 




:e = 


11.206850 
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APPENDIX F 



SAMPLE INPUT/OUTPUT DATA FILES USED WITH NEC 



This Appendix lists an input data file used NEC and the 
corresponding output file. The run selected was for an 
antenna with the following physical characteristics using one 
reflection: 



Driving Frequency 
Antenna Height 
Antenna Radius 
Gap Distance 
Ground Plate Spacing 



299.8 MHz 
0.48 wavelengths 
0.01 wavelengths 
0.02 wavelenghts 
1.4 wavelengths 



CE 




Run 


11 


Monopole Antenna of ! 


GW 


i, 


24, 


0, 


0, 0, 0, 


0, 


.48, .01 


GW 


2, 


48, 


0, 


0, 2.32, 


0, 


0, 3.28, . 


GE 


1 












GN 


1 












EK 














EX 


0 


1 




1 0 




1.0 


EX 


0 


2 




24 0 




-1.0 


EX 


0 


2 




25 0 




-1.0 


XQ 














EN 
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Sample output file from the above data is shown below: 



************************************ 

NUMERICAL ELECTROMAGNETICS CODE 
************************************ 



COMMENTS 



Run 11: Monopole Antenna of Height .48, Radius .01,. 02, 1.40 



STRUCTURE SPECIFICATION 

COORDINATES MUST BE INPUT IN 
METERS OR BE SCALED TO METERS 
BEFORE STRUCTURE INPUT IS ENDED 



WIRE 

LAST 


TAG 














NO. OF 


FIRST 


NO. 

SEG. 


XI 

NO. 


Y1 


Z1 


X2 


Y2 


Z2 


RADIUS 


SEG. 


SEG. 


1 

24 


.00000 

1 


.00000 


.00000 


.00000 


.00000 


.48000 


.01000 


24 


1 


2 

72 


.00000 

2 


.00000 


2.32000 


.00000 


.00000 


3.28000 


.01000 


48 


25 



GROUND PLANE SPECIFIED. 

WHERE WIRE ENDS TOUCH GROUND, CURRENT WILL BE INTERPOLATED TO IMAGE IN GROUND PLANE. 



TOTAL SEGMENTS USED- 72 NO. SEG. IN A SYMMETRIC CELL* 72 SYMMETRY FLAG= 0 



- MULTIPLE WIRE JUNCTIONS - 
JUNCTION SEGMENTS (- FOR END 1, + FOR END 2) 
NONE 



SEGMENTATION DATA 

COORDINATES IN METERS 

1+ AND I- INDICATE THE SEGMENTSBEFORE AND AFTER I 



SEG. 


COORDINATES 


OF SEG. 


CENTER 


SEG. 


ORIENTATION 


ANGLES 


WIRE 


CONNECTION 


DATA 


TAG 


NO. 


X 


Y 


z 


LENGTH 


ALPHA 


BETA 


RADIUS 


I- 


I 


1+ 


NO. 


1 


.00000 


.00000 


.01000 


.02000 


90.00000 


.00000 


.01000 


1 


1 


2 


1 


2 


.00000 


.00000 


.03000 


.02000 


90.00000 


.00000 


.01000 


1 


2 


3 


1 


3 


.00000 


.00000 


.05000 


.02000 


90.00000 


.00000 


.01000 


2 


3 


4 


1 


4 


.00000 


.00000 


.07000 


.02000 


90.00000 


.00000 


.01000 


3 


4 


5 


1 


5 


.00000 


.00000 


.09000 


.02000 


90.00000 


.00000 


.01000 


4 


5 


6 


1 


6 


.00000 


.00000 


.11000 


.02000 


90.00000 


.00000 


.01000 


5 


6 


7 


1 


7 


.00000 


.00000 


.13000 


.02000 


90.00000 


.00000 


.01000 


6 


7 


8 


1 


8 


.00000 


.00000 


.15000 


.02000 


90.00000 


.00000 


.01000 


7 


8 


9 


1 


9 


.00000 


.00000 


.17000 


.02000 


90.00000 


.00000 


.01000 


8 


9 


10 


1 


10 


.00000 


.00000 


.19000 


.02000 


90.00000 


.00000 


.01000 


9 


10 


11 


1 


11 


.00000 


.00000 


.21000 


.02000 


90.00000 


.00000 


.01000 


10 


11 


12 


1 


12 


.00000 


.00000 


.23000 


.02000 


90.00000 


.00000 


.01000 


11 


12 


13 


1 


13 


.00000 


.00000 


.25000 


.02000 


90.00000 


.00000 


.01000 


12 


13 


14 


1 



109 



14 


.00000 


.00000 


.27000 


.02000 


90.00000 


.00000 


.01000 


13 


14 


15 


1 


15 


.00000 


.00000 


.29000 


.02000 


90.00000 


.00000 


.01000 


14 


15 


16 


1 


16 


.00000 


.00000 


.31000 


.02000 


90.00000 


.00000 


.01000 


15 


16 


17 


1 


17 


.00000 


.00000 


.33000 


.02000 


90.00000 


.00000 


.01000 


16 


17 


18 


1 


18 


.00000 


.00000 


.35000 


.02000 


90.00000 


.00000 


.01000 


17 


18 


19 


1 


19 


.00000 


.00000 


.37000 


.02000 


90.00000 


.00000 


.01000 


18 


19 


20 


1 


20 


.00000 


.00000 


.39000 


.02000 


90.00000 


.00000 


.01000 


19 


20 


21 


1 


21 


.00000 


.00000 


.41000 


.02000 


90.00000 


.00000 


.01000 


20 


21 


22 


1 


22 


.00000 


.00000 


.43000 


.02000 


90.00000 


.00000 


,01000 


21 


22 


23 


1 


23 


.00000 


.00000 


.45000 


.02000 


90.00000 


.00000 


,01000 


22 


23 


24 


1 


24 


.00000 


.00000 


.47000 


.02000 


90.00000 


.00000 


.01000 


23 


24 


0 


1 


25 


.00000 


.00000 


2.33000 


.02000 


90.00000 


.00000 


.01000 


0 


25 


26 


2 


26 


.00000 


.00000 


2.35000 


.02000 


90.00000 


.00000 


.01000 


25 


26 


27 


2 


27 


.00000 


.00000 


2.37000 


.02000 


90.00000 


.00000 


.01000 


26 


27 


28 


2 


28 


.00000 


.00000 


2.39000 


.02000 


90.00000 


.00000 


.01000 


27 


28 


29 


2 


29 


.00000 


.00000 


2.41000 


.02000 


90.00000 


.00000 


.01000 


28 


29 


30 


2 


30 


.00000 


.00000 


2.43000 


.02000 


90.00000 


.00000 


.01000 


29 


30 


31 


2 


31 


.00000 


.00000 


2.45000 


.02000 


90.00000 


.00000 


.01000 


30 


31 


32 


2 


32 


.00000 


.00000 


2.47000 


.02000 


90.00000 


.00000 


.01000 


31 


32 


33 


2 


33 


.00000 


.00000 


2.49000 


.02000 


90.00000 


.00000 


.01000 


32 


33 


34 


2 


34 


.00000 


.00000 


2.51000 


.02000 


90.00000 


.00000 


.01000 


33 


34 


35 


2 


35 


.00000 


.00000 


2.53000 


.02000 


90.00000 


.00000 


.01000 


34 


35 


36 


2 


36 


.00000 


.00000 


2.55000 


.02000 


90.00000 


.00000 


.01000 


35 


36 


37 


2 


37 


.00000 


.00000 


2.57000 


.02000 


90.00000 


.00000 


.01000 


36 


37 


38 


2 


38 


.00000 


.00000 


2.59000 


.02000 


90.00000 


.00000 


.01000 


37 


38 


39 


2 


39 


.00000 


.00000 


2.61000 


.02000 


90.00000 


.00000 


.01000 


38 


39 


40 


2 


40 


.00000 


.00000 


2.63000 


.02000 


90.00000 


.00000 


.01000 


39 


40 


41 


2 


41 


.00000 


.00000 


2.65000 


.02000 


90.00000 


.00000 


.01000 


40 


41 


42 


2 


42 


.00000 


.00000 


2.67000 


.02000 


90.00000 


.00000 


.01000 


41 


42 


43 


2 


43 


.00000 


.00000 


2.69000 


.02000 


90.00000 


.00000 


.01000 


42 


43 


44 


2 


44 


.00000 


.00000 


2.71000 


.02000 


90.00000 


.00000 


.01000 


43 


44 


45 


2 


45 


.00000 


.00000 


2.73000 


.02000 


90.00000 


.00000 


.01000 


44 


45 


46 


2 


46 


.00000 


.00000 


2.75000 


.02000 


90.00000 


.00000 


.01000 


45 


46 


47 


2 


47 


.00000 


.00000 


2.77000 


.02000 


90.00000 


.00000 


.01000 


46 


47 


48 


2 


48 


.00000 


.00000 


2.79000 


.02000 


90.00000 


.00000 


.01000 


47 


48 


49 


2 


49 


.00000 


.00000 


2.81000 


.02000 


90.00000 


.00000 


.01000 


48 


49 


50 


2 


50 


.00000 


.00000 


2.83000 


.02000 


90.00000 


.00000 


.01000 


49 


50 


51 


2 


51 


.00000 


.00000 


2.85000 


.02000 


90.00000 


.00000 


.01000 


50 


51 


52 


2 


52 


.00000 


.00000 


2.87000 


.02000 


90.00000 


.00000 


.01000 


51 


52 


53 


2 


53 


.00000 


.00000 


2.89000 


.02000 


90.00000 


.00000 


.01000 


52 


53 


54 


2 


54 


.00000 


.00000 


2.91000 


.02000 


90.00000 


.00000 


.01000 


53 


54 


55 


2 


55 


.00000 


.00000 


2.93000 


.02000 


90.00000 


.00000 


.01000 


54 


55 


56 


2 


56 


.00000 


.00000 


2.95000 


.02000 


90.00000 


.00000 


.01000 


55 


56 


57 


2 


57 


.00000 


.00000 


2.97000 


.02000 


90.00000 


.00000 


.01000 


56 


57 


58 


2 


58 


.00000 


.00000 


2.99000 


.02000 


90.00000 


.00000 


.01000 


57 


58 


59 


2 


59 


.00000 


.00000 


3.01000 


.02000 


90.00000 


.00000 


.01000 


58 


59 


60 


2 


60 


.00000 


.00000 


3.03000 


.02000 


90.00000 


.00000 


.01000 


59 


60 


61 


2 


61 


.00000 


.00000 


3.05000 


.02000 


90.00000 


.00000 


.01000 


60 


61 


62 


2 


62 


.00000 


.00000 


3.07000 


.02000 


90.00000 


.00000 


.01000 


61 


62 


63 


2 


63 


.00000 


.00000 


3.09000 


.02000 


90.00000 


.00000 


.01000 


62 


63 


64 


2 


64 


.00000 


.00000 


3.11000 


.02000 


90.00000 


.00000 


.01000 


63 


64 


65 


2 


65 


.00000 


.00000 


3.13000 


.02000 


90.00000 


.00000 


.01000 


64 


65 


66 


2 


66 


.00000 


.00000 


3.15000 


.02000 


90.00000 


.00000 


.01000 


65 


66 


67 


2 


67 


.00000 


.00000 


3.17000 


.02000 


90.00000 


.00000 


.01000 


66 


67 


68 


2 


68 


.00000 


.00000 


3.19000 


.02000 


90.00000 


.00000 


.01000 


67 


68 


69 


2 


69 


.00000 


.00000 


3.21000 


.02000 


90.00000 


.00000 


.01000 


68 


69 


70 


2 


70 


.00000 


.00000 


3.23000 


.02000 


90.00000 


.00000 


.01000 


69 


70 


71 


2 


71 


.00000 


.00000 


3.25000 


.02000 


90.00000 


.00000 


.01000 


70 


71 


72 


2 


72 


.00000 


.00000 


3.27000 


.02000 


90.00000 


.00000 


.01000 


71 


72 


0 


2 



***** DATA 
0.00000 E +00 


CARD NO . 1 
0 . 00000 E +00 


GN 


1 


0 


0 


0 0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


***** DATA 
0 . 00000 E +00 


CARD NO . 2 
0 . OOOOOE+OO 


EK 


0 


0 


0 


0 0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


***** DATA 
0 . 00000 E +00 


CARD NO . 3 
0. 00000 E +00 


EX 


0 


1 


1 


0 1. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


***** DATA 
0.00000 E +00 


CARD NO . 4 
0. OOOOOE+OO 


EX 


0 


2 


24 


0 -1. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


***** DATA 
0 . 00000 E +00 


CARD NO . 5 
0. OOOOOE+OO 


EX 


0 


2 


25 


0 -1. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 


0. OOOOOE+OO 
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***** DATA CARD NO. 6 XQ 0 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 

0 . 00000E+00 0 . 00000E+00 



• FREQUENCY 

FREQUENCY- 2.9980E+02 MHZ 
WAVELENGTH- 1.0000E+00 METERS 



APPROXIMATE INTEGRATION EMPLOYED FOR SEGMENTS MORE THAN 1.000 WAVELENGTHS APART 
THE EXTENDED THIN WIRE KERNEL WILL BE USED 

STRUCTURE IMPEDANCE LOADING 

THIS STRUCTURE IS NOT LOADED 

ANTENNA ENVIRONMENT 

PERFECT GROUND 



MATRIX TIMING 

FILL- 1.500 MIN., FACTOR- .364 MIN. 



ANTENNA INPUT PARAMETERS 

TAG SEG. VOLTAGE (VOLTS) CURRENT (AMPS) IMPEDANCE (OHMS) ADMITTANCE 

(MHOS) POWER 

NO. NO. REAL IMAG. REAL IMAG. REAL IMAG. REAL 

IMAG. (WATTS) 

1 1 1 . 00000E+00 0 . 00000E+00 2.43054E-03 2.78753E-03 1 . 77699E+02-2 . 03799E+02 2.43054E-03 
2.78753E-03 1.21527E-03 

2 48-1.00000E+00 0 . OOOOOE+OO-2 . 42866E-03-2. 79143E-03 1 . 77398E+02-2 . 03896E+02 2.42866E-03 
2.79143E-03 1.21433E-03 

2 49-1.00000E+00 0 . 00000E+00-2. 42265E-03-2. 79230E-03 1 . 77274E+02-2 . 04322E+02 2.42265E-03 

2.79230E-03 1.21132E-03 



CURRENTS AND LOCATION 

DISTANCES IN WAVELENGTHS 



SEG. 


TAG 


COORD. 


OF SEG. 


CENTER 


SEG. 


CURRENT (AMPS) 




NO. 


NO. 


X 


Y 


z 


LENGTH 


REAL IMAG. 


MAG. 


PHASE 


1 


1 


.0000 


.0000 


.0100 


.02000 


2.4305E-03 2.7875E-03 


3.6984E-03 


48.914 


2 


1 


.0000 


.0000 


.0300 


.02000 


2.4147E-03 1.0516E-03 


2.6338E-03 


23.534 


3 


1 


.0000 


.0000 


.0500 


.02000 


2.3833E-03 -8.5797E-05 


2.3848E-03 


-2.062 


4 


1 


.0000 


.0000 


.0700 


.02000 


2.3366E-03 -9.8162E-04 


2.5344E-03 


-22.787 


5 


1 


.0000 


.0000 


.0900 


.02000 


2.2754E-03 -1.7685E-03 


2.8819E-03 


-37.856 


6 


1 


.0000 


.0000 


.1100 


.02000 


2.2004E-03 -2.4668E-03 


3.3055E-03 


-48.267 


7 


1 


.0000 


.0000 


.1300 


.02000 


2.1125E-03 -3.0847E-03 


3 .7387E-03 


-55.595 


8 


1 


.0000 


.0000 


.1500 


.02000 


2.0130E-03 -3.6240E-03 


4 . 1456E-03 


-60.950 


9 


1 


.0000 


.0000 


.1700 


.02000 


1.9031E-03 -4.0841E-03 


4.5058E-03 


-65.016 


10 


1 


.0000 


.0000 


.1900 


.02000 


1.7842E-03 -4.4630E-03 


4.8065E-03 


-68.209 


11 


1 


.0000 


.0000 


.2100 


.02000 


1.6579E-03 -4.7589E-03 


5.0394E-03 


-70.792 


12 


1 


.0000 


.0000 


.2300 


.02000 


1.5258E-03 -4.9697E-03 


5. 1987E-03 


-72.932 


13 


1 


.0000 


.0000 


.2500 


.02000 


1.3895E-03 -5.0945E-03 


5.2806E-03 


-74.744 


14 


1 


.0000 


.0000 


.2700 


.02000 


1.2508E-03 -5.1325E-03 


5.2827E-03 


-76.304 


15 


1 


.0000 


.0000 


.2900 


.02000 


1.1113E-03 -5.0841E-03 


5.2041E-03 


-77.670 


16 


1 


.0000 


.0000 


.3100 


.02000 


9.7278E-04 -4.9502E-03 


5.0449E-03 


-78.882 
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17 


1 


.0000 


.0000 


.3300 


.02000 


8.3692E-04 


-4.7327E-03 


4.8061E-03 


-79.972 


18 


1 


.0000 


.0000 


.3500 


.02000 


7.0529E-04 


-4.4340E-03 


4.4897E-03 


-80.962 


19 


1 


.0000 


.0000 


.3700 


.02000 


5.7941E-04 


-4.0569E-03 


4.0981E-03 


-81.872 


20 


1 


.0000 


.0000 


.3900 


.02000 


4.6064E-04 


-3.6044E-03 


3 . 6337E-03 


-82.717 


21 


1 


.0000 


.0000 


.4100 


.02000 


3.5018E-04 


-3.0785E-03 


3.0983E-03 


-83.510 


22 


1 


.0000 


.0000 


.4300 


.02000 


2.4894E-04 


-2.4785E-03 


2.4909E-03 


-84.264 


23 


1 


.0000 


.0000 


.4500 


.02000 


1.5733E-04 


-1.7963E-03 


1.8032E-03 


-84.995 


24 


1 


.0000 


.0000 


.4700 


.02000 


7.2003E-05 


-9.5821E-04 


9.6091E-04 


-85.703 


25 


2 


.0000 


.0000 


2.3300 


.02000 


-8.2718E-05 


9.5969E-04 


9.6325E-04 


94.926 


26 


2 


.0000 


.0000 


2.3500 


.02000 


-1.7730E-04 


1.7990E-03 


1.8077E-03 


95.629 


27 


2 


.0000 


.0000 


2.3700 


.02000 


-2.7635E-04 


2.4819E-03 


2.4972E-03 


96.353 


28 


2 


.0000 


.0000 


2.3900 


.02000 


-3.8402E-04 


3.0826E-03 


3.1064E-03 


97.101 


29 


2 


.0000 


.0000 


2.4100 


.02000 


-5.0002E-04 


3.6090E-03 


3.6434E-03 


97.888 


30 


2 


.0000 


.0000 


2.4300 


.02000 


-6.2347E-04 


4.0618E-03 


4.1094E-03 


98.727 


31 


2 


.0000 


.0000 


2.4500 


.02000 


-7.5316E-04 


4.4391E-03 


4 . 5025E-03 


99.629 


32 


2 


.0000 


.0000 


2.4700 


.02000 


-8.8769E-04 


4.7379E-03 


4.8203E-03 


100.612 


33 


2 


.0000 


.0000 


2.4900 


.02000 


-1.0256E-03 


4.9554E-03 


5.0604E-03 


101.693 


34 


2 


.0000 


.0000 


2.5100 


.02000 


-1.1652E-03 


5.0891E-03 


5. 2208E-03 


102.896 


35 


2 


.0000 


.0000 


2.5300 


.02000 


-1.3049E-03 


5.1373E-03 


5. 3004E-03 


104.252 


36 


2 


.0000 


.0000 


2.5500 


.02000 


-1.4430E-03 


5. 0989E-03 


5.2992E-03 


105.801 


37 


2 


.0000 


.0000 


2.5700 


.02000 


-1.5778E-03 


4.9737E-03 


5.2180E-03 


107.600 


38 


2 


.0000 


.0000 


2.5900 


.02000 


-1.7076E-03 


4.7623E-03 


5.0592E-03 


109.726 


39 


2 


.0000 


.0000 


2.6100 


.02000 


-1.8308E-03 


4.4659E-03 


4.8266E-03 


112.292 


40 


2 


.0000 


.0000 


2.6300 


.02000 


-1.9460E-03 


4.0864E-03 


4.5261E-03 


115.464 


41 


2 


.0000 


.0000 


2.6500 


.02000 


-2.0515E-03 


3.6256E-03 


4.1658E-03 


119.502 


42 


2 


.0000 


.0000 


2.6700 


.02000 


-2.1461E-03 


3.0856E-03 


3.7585E-03 


124.820 


43 


2 


.0000 


.0000 


2.6900 


.02000 


-2.2286E-03 


2.4669E-03 


3.3245E-03 


132.094 


44 


2 


.0000 


.0000 


2.7100 


.02000 


-2. 2979E-03 


1 . 7679E-03 


2.8993E-03 


142.427 


45 


2 


.0000 


.0000 


2.7300 


.02000 


-2.3532E-03 


9.8018E-04 


2.5492E-03 


157.387 


46 


2 


.0000 


.0000 


2.7500 


.02000 


-2.3937E-03 


8.3523E-05 


2.3952E-03 


178.002 


47 


2 


.0000 


.0000 


2.7700 


.02000 


-2.4190E-03 


-1.0547E-03 


2.6389E-03 


-156.441 


48 


2 


.0000 


.0000 


2.7900 


.02000 


-2.4287E-03 


-2.7914E-03 


3.7001E-03 


-131.025 


49 


2 


.0000 


.0000 


2.8100 


.02000 


-2.4226E-03 


-2.7923E-03 


3.6968E-03 


-130.945 


50 


2 


.0000 


.0000 


2.8300 


.02000 


-2.4010E-03 


-1.0573E-03 


2.6235E-03 


-156.235 


51 


2 


.0000 


.0000 


2.8500 


.02000 


-2.3641E-03 


7.9330E-05 


2.3655E-03 


178.078 


52 


2 


.0000 


.0000 


2.8700 


.02000 


-2.3124E-03 


9.7436E-04 


2.5093E-03 


157.151 


53 


2 


.0000 


.0000 


2.8900 


.02000 


-2.2466E-03 


1.7605E-03 


2.8542E-03 


141.917 


54 


2 


.0000 


.0000 


2.9100 


.02000 


-2 . 1675E-03 


2.4580E-03 


3.2772E-03 


131.406 


55 


2 


.0000 


.0000 


2.9300 


.02000 


-2.0762E-03 


3.0752E-03 


3.7104E-03 


124.025 


56 


2 


.0000 


.0000 


2.9500 


.02000 


-1.9738E-03 


3.6139E-03 


4 .1178E-03 


118.642 


57 


2 


.0000 


.0000 


2.9700 


.02000 


-1.8617E-03 


4.0734E-03 


4.4786E-03 


114.562 


58 


2 


.0000 


.0000 


2.9900 


.02000 


-1.7413E-03 


4.4518E-03 


4.7803E-03 


111.362 


59 


2 


.0000 


.0000 


3.0100 


.02000 


-1.6141E-03 


4.7472E-03 


5.0141E-03 


108.779 


60 


2 


.0000 


.0000 


3.0300 


.02000 


-1.4819E-03 


4 . 9578E-03 


5.1745E-03 


106.641 


61 


2 


.0000 


.0000 


3.0500 


.02000 


-1.3461E-03 


5.0823E-03 


5.2575E-03 


104.834 


62 


2 


.0000 


.0000 


3.0700 


.02000 


-1.2084E-03 


5.1202E-03 


5.2609E-03 


103.280 


63 


2 


.0000 


.0000 


3.0900 


.02000 


- 1 . 0707E-03 


5. 0718E-03 


5.1836E-03 


101.921 


64 


2 


.0000 


.0000 


3.1100 


.02000 


-9.3452E-04 


4.9381E-03 


5.0257E-03 


100.716 


65 


2 


.0000 


.0000 


3.1300 


.02000 


-8.0147E-04 


4.7209E-03 


4.7884E-03 


99.635 


66 


2 


.0000 


.0000 


3.1500 


.02000 


-6.7310E-04 


4.4227E-03 


4.4736E-03 


98.654 


67 


2 


.0000 


.0000 


3.1700 


.02000 


-5.5087E-04 


4.0463E-03 


4.0836E-03 


97.753 


68 


2 


.0000 


.0000 


3.1900 


.02000 


-4.3609E-04 


3.5946E-03 


3.6210E-03 


96.917 


69 


2 


.0000 


.0000 


3.2100 


.02000 


-3.2990E-04 


3.0698E-03 


3.0875E-03 


96.134 


70 


2 


.0000 


.0000 


3.2300 


.02000 


-2.3317E-04 


2.4712E-03 


2.4822E-03 


95.390 


71 


2 


.0000 


.0000 


3.2500 


.02000 


-1.4630E-04 


1.7909E-03 


1.7969E-03 


94.670 


72 


2 


.0000 


.0000 


3.2700 


.02000 


-6.6339E-05 


9.5517E-04 


9 . 5747E-04 


93.973 



POWER BUDGET 

INPUT POWER - 3.6409E-03 WATTS 
RADIATED POWER- 3.6409E-03 WATTS 
STRUCTURE LOSS- 0.0000E+00 WATTS 
NETWORK LOSS - 0.0000E+00 WATTS 
EFFICIENCY - 100.00 PERCENT 



***** DATA CARD NO. 7 EN 0 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 

0 . 00000E+00 0 . 00000E+00 

RUN TIME = 1.952 
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